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00 Abstract 

r~i We review the recent developments in applying holographic methods to understand 

non-equilibrium physics in strongly coupled field theories. The emphasis will be on 
^ elucidating the relation between evolution of quantum field theories perturbed away 

4^3 from equilibrium and the dual picture of dynamics of classical fields in black hole 

backgrounds. In particular, we discuss the linear response regime, the hydrodynamic 
regime and finally the non- linear regime of interacting quantum systems. We also 
describe how the duality might be used to learn some salient aspects of black hole 
physics in terms of field theory observables. 
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1 Introduction 

Astounding amount of progress and understanding in physics has been achieved by studying 
special systems in equilibrium, which are by definition 'non-dynamical', i.e. independent of 
time. One important reason for the prevalent focus on time-independent systems is that 
they are of course far simpler to study than dynamical ones. Indeed, vast majority of ex- 
actly tractable systems are of this kind. Known exact solutions are seldom fully generic 
and typically admit a large degree of symmetry, usually including time-translation invari- 
ance. Nevertheless, while the study of such non-dynamical, stationary situations might seem 
rather looking-under-the-lamppost type strategy, it has proved to be a very successful one. 
Although due to an incredible complexity of nature around us such exactly-tractable systems 
are at best only approximations to the real world, they are often remarkably relevant and 
useful. Dynamically-evolving systems tend to equilibrate, and in absence of external forcing 
they typically settle down to stationary configurations.^ Thus stable stationary solutions can 
reveal the late-time physics of a generic system. Moreover, many interesting physical pro- 
cesses such as phase transitions, which we observe occurring dynamically, can be well-studied 
without invoking any explicit time-dependence at all. 

Nevertheless, not many would dispute that most interesting physical phenomena do in- 
volve non-trivial temporal dynamics. Not only are dynamically evolving systems ubiquitous 
in nature, but they are the raison d'etre for everything we see around us. Hence the need 
to understand dynamics scarcely needs motivation. However, progress so far has unfortu- 
nately been hindered by lack of adequate techniques. Typically, one resorts to perturbative 

^ Implicit in this statement is an assumption that the systems under consideration are sufficiently ergodic. 
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methods, though the regime of their vahdity places severe hmitations on their apphcabihty. 
Alternately, one might try to "put the system on a computer" and evolve numerically, but 
the computational cost involved is usually too astronomical to allow for convenient extraction 
of the physics. In certain cases one may circumvent both of these limitations by mapping 
the system into a more tractable one. We will see that all these ingredients come into play 
in the present Review. 

So far our somewhat self-evident remarks have been rather abstract and general. We 
will now specify the particular context we wish to address. We will concentrate on exploring 
the dynamics of a certain class of quantum field theories, focusing especially on the strong 
coupling regime. Although we have best understanding of field theories at weak coupling, 
strongly-coupled systems are ubiquitous in nature, ranging from typical condensed matter 
systems, to quark-gluon plasmas created in high-energy experiments, and in fact play a role 
in most areas of physics. Apart from the evident applicability, there is also the rewarding 
aspect of serendipity related to the present cutting-edge experiments involving these systems. 

Which properties of strongly coupled field theories do we wish to understand, and which 
ones can we hope to understand? Fortuitously there is a substantial overlap between the 
answers to both questions. In particular, it is both interesting and tractable to extract 
certain universal features, as we will discuss below. Conversely, we seldom can, nor want, 
to calculate the detailed microscopic dynamics, due to the sheer level of complexity; it is 
generally much more instructive to take a coarse-grained view of the system. As a result we 
will often focus on obtaining a low energy effective description for such strongly interacting 
systems. From a conventional renormalization group (RG) picture, it is then clear that one 
expects similar low energy physics for systems within the same universality class. One such 
low energy theory which we will discuss in some detail is hydrodynamics, or more generally 
fluid dynamics, which is expected to be a good description as long as the local fluid variables 
vary slowly compared to the microscopic scale, i.e. at long wavelengths and small frequencies. 

Hence in studying dynamics of strongly coupled field theories, we are simultaneously 
exploring the dynamics of fluids. Our study is then bolstered by the insights which we have 
already acquired from hydrodynamics. On the other hand, despite decades of theoretical 
as well as numerical, observational, and experimental scrutiny which fluid dynamics has re- 
ceived, there are still many deep questions which remain to be answered, especially involving 
dynamical evolution. For example, one of the famous Clay Millennium Prize Problems con- 
cerns the global regularity (existence and smoothness) of the Navier-Stokes equations [1]. 
Intriguingly, the solutions often include turbulence, which, in spite of its practical importance 
in science and engineering, still remains one of the great unsolved problems in physics. 

As already mentioned above, understanding dynamics in strongly coupled field theories, 
or their effective description in terms of fluids, is an exceedingly hard problem. One of 
the key strategies has been to focus on field theories which admit a holographic dual, and 
use this dual description to extract the physical properties of the field theoretic system 
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under consideration. The prototypical case is the AdS/CFT correspondence [2, 3, 4] which 
relates the four-dimensional A/" = 4 Super Yang-Mills (SYM) gauge theory to a IIB string 
theory (or supergravity) on asymptotically AdSs x spacetime.^ As is well known, this 
is a strong/weak coupling duality; the strongly-coupled field theory can be accessed via the 
semi-classical gravitational dual, which has obvious computational, as well as conceptual, 
advantages. If we wish to know how a given strongly-coupled system behaves, we need 
not devise methods to calculate this in the field theory at all - a rather daunting, if not 
impossible, task - we only need to translate the system into the dual language and calculate 
the classical evolution using Einstein's equations. 

While the Af = 4 SYM is fundamentally distinct from the 'real-world' systems which we 
would ultimately like to understand, it can serve as a useful toy model to motivate and study 
new classes of strongly coupled phenomena. This philosophy has led to the correspondence 
being applied to QCD via the so called AdS/QCD approach (for reviews see [8, 9, 10, 11]), 
and more recently to condensed matter systems, often dubbed AdS/CMT, where a variety 
of physical effects raging from superfluid transitions to non- Fermi liquid behaviour are being 
actively investigated. An excellent account of these efforts can be found in the reviews 
[12, 13, 7, 14]. Of course, what makes this enterprise fascinating is the degree to which the 
computations in the AdS/CFT framework agree with experimentally measured quantities in 
the real world. One such quantity which received much attention is the shear viscosity of the 
quark-gluon plasma; based holographic computations, [15] suggest that the dimensionless 
ratio of viscosity to entropy density {i]/s) has a universal lower bound, i]/s > l/47r, which 
subsequently has been subject to intense scrutiny; for the current status in string theory see 
[16, 17]. 

Perhaps the most intriguing - and promising - aspect of the gauge/gravity correspondence 
is highlighted by its holographic nature: the theories which are dual to each other are 
naturally formulated in different number of dimensions.^ This automatically implies that 
the effective degrees of freedom on the two sides of the correspondence are related in a 
highly non-local (and hitherto rather mysterious) fashion. This of course makes the task 
of extracting physics of one theory from its dual formulation even more formidable than 
merely doing calculations in the dual; indeed, much of the AdS/CFT-related efforts of the 
last decade have concentrated on elucidating the dictionary between the two sides. Yet at the 
same time this very feature, which might initially appear as an unwanted complication, lies 
at the heart of the enormous potential the duality holds for solving the system. The effects 
we seek to unravel are often very complex, emergent phenomena which cannot (in practice, 
or even in principle) be described in terms of the fundamental degrees of freedom. In other 
words, they require a different description. This by itself does not of course guarantee that 

^ The AdS/CFT correspondence is comprehensively reviewed in the classic reviews [5, 6]. See also [7] for 
a nice review with emphasis on 'applied holography'. 

^ The holographic principle was proposed much before the advent of the AdS/CFT correspondence and 
was motivated in part by the peculiar non-extensive nature of black hole entropy in [18, 19]. 
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the holographic description is the desired one, but the fact that out of such formidably 
complicated repackaging of the information a new simple and elegant picture arises, is at 
least promising. It may well transpire that the complex emergent phenomena we seek to 
understand, at least those which are in some sense fundamental, appear simple when re- 
written in the dual language. 

One classic indication that the above hope is more that just a wishful thinking may be 
found in the manner in which the radial direction of the bulk emerges from the boundary field 
theory, sometimes referred to as the scale/radius, or UV/IR, duality [20]. From everyday 
experience, we are well-aware that physics likes to organize itself by energy (or length) scales; 
processes occurring at widely-separated scales do not interact very much. In the holographic 
dual, this hierarchy of scales, mathematically expressed in terms of Wilsonian RG, is neatly 
packaged in terms of an emergent radial direction of the bulk geometry (this idea is at the 
heart of the holographic renormalization group developed initially in [21]). But rather than 
just serving as a mnemonic for the energy scale, this dimension takes on a life of its own: it 
mixes with the other directions in a fully covariant fashion! Taking the bulk perspective, the 
field theory's decoupling of scales is simply a manifestation of bulk locality. So this simple 
and natural feature of the bulk description has far-reaching implications for the boundary 
description. 

Phrased more prosaically, it is not just to be hoped-for, but rather it is guaranteed, 
that certain complex phenomena must have a very simple and natural explanation in the 
dual picture. After all, the two dual theories are just different descriptions of the same 
physical system: there is no absolute notion of which side is 'more fundamental' than the 
other. Whichever side we use as a starting point, there are certain fundamental properties 
or principles underlying the theory that we understand in simple terms; yet these will be 
mapped into the dual framework, wherein they will appear as highly non-trivial. If we 
were to come upon them from the other side - which due to their inherent importance is 
perhaps not entirely unlikely - it would appear rather magical that in the dual language they 
suddenly become simple. Of course, this is not to say that all complex phenomena must 
admit a simple description in some reformulation, but it is not so outrageous to expect that 
the fundamentally important ones do. 

The preceding philosophical interlude was meant to motivate the use of gauge/gravity 
duality to understand certain field theories at strong coupling in terms of their gravitational 
dual, beyond the mere fact that we can usually calculate in classical gravity more easily 
than in the strongly-coupled field theory. But the gauge/gravity correspondence likewise 
has profound implications when applied in the other direction. The field theory, albeit 
strongly coupled, provides a definition of quantum gravity with asymptotically AdS boundary 
conditions. Hence even in extreme regimes where classical gravity breaks down, and where 
we don't yet have the tools to understand the physics within string theory directly, the field 
theory remains well-defined. This means that once we achieve sufficient understanding of the 
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dictionary between the two sides, we can elucidate the long-standing quantum gravitational 
puzzles by re-casting them into the field theoretic language. 

Let us briefiy indicate some of the gravitational questions we ultimately hope to answer, 
since these have posed the underlying motivation in many of the works mentioned in this 
Review. A central question of quantum gravity concerns the fundamental nature of space- 
time. We have come to realize that spacetime is an emergent concept, but exactly how it 
emerges remains a mystery. Happily, certain aspects of this emergence can be conveniently 
explored using the AdS/CFT framework, including ones pertaining to the present theme of 
out-of-equilibrium dynamics. 

To set the stage, let us first recall several well-understood classical highlights of the corre- 
spondence. According to the AdS/CFT dictionary, different asymptotically AdS spacetimes 
manifest themselves by different states in the boundary field theory. For example, the vac- 
uum state in the CFT corresponds to the pure AdS spacetime. Metric perturbations which 
maintain the AdS asymtopia are related to the stress-energy-momentum tensor expectation 
value in the CFT. More importantly, putting a black hole in the bulk has the effect of heating 
up the boundary theory. Specifically, a large^ Schwarzschild-AdS black hole corresponds to 
(approximately) thermal state in the gauge theory. This can be easily conceptualized as the 
late-time configuration a generic state evolves to: in the bulk, the combined effect of gravity 
and negative curvature tends to make a generic large-energy configuration collapse and form 
a black hole which then quickly settles down to the Schwarzschild-AdS geometry. On the 
other hand, in the field theory, a generic large-energy excitation will thermalize. At the level 
of this coarse entry in the dictionary we see that heating the field theory corresponds to 
black hole formation in the dual gravity, and the subsequent thermalization corresponds to 
the black hole settling down to a stationary state. 

While appealingly simple, this level of understanding is far too coarse to allow us to 
extract the more interesting aspects. We need to probe the AdS/CFT dictionary further 
to uncover what happens in regions where the classical description of the black hole breaks 
down, such as near the curvature singularity, or in the more general dynamical situations. 
Ultimately, we would like to answer such questions as: Which CFT configurations admit a 
dual description in terms of classical spacetime? What types of spacetime singularities are 
physically allowed? How are the disallowed singularities resolved? How is spacetime causal 
structure encoded in the dual field theory? 

Emergence of time is, if anything, even more mysterious than the emergence of space. 
Not only do we have more satisfactory toy models of the latter than of the former, but con- 
ceptually the problem of time is one of the deepest problems in quantum gravity. The 'time' 
(conjugate to the Hamiltonian) which quantum mechanics uses for evolution is ingrained 

^ AdS is a space of constant negative curvature, which introduces a length scale, called the AdS scale 
-RAdSj corresponding to the radius of curvature. The black hole size is then measured in terms of this AdS 
scale; large black holes have horizon radius r+ > i?Ads- 
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in the fundamental formulation of the system; it is there from the start rather than being 
emergent. We can take this evolution parameter to be simply the time of the non-dynamical 
(and non-emergent) background on which the CFT lives. What is then the relation between 
this CFT time and the notion of time in the dual bulk spacetime? For bulk spacetimes which 
are globally static we tend to associate the two; but already for general dynamical space- 
times there is no natural identification even classically, since there is no uniquely-specified 
foliation of the bulk. Even if we can construct geometrically defined spacelike slices through 
the bulk (such as volume- maximizing ones), there is a-priori no reason that the bulk events 
localized on a given slice should be dual to boundary events localized at the correspond- 
ing boundary time; to the contrary, we have many indications that the correspondence is 
much more temporally non-local. Nevertheless, we still expect that time-dependence on the 
boundary will be manifested by time-dependence in the bulk, and vice-versa. Indeed, focus- 
ing on certain characteristic features of a given time-dependence can enable us to elucidate 
the gauge/gravity dictionary further. Thus, our overarching motif of considering dynamical 
systems in strongly coupled field theories will naturally translate to studying bulk spacetime 
dynamics. 

We now return to our initial comments regarding time-dependence being difficult to 
handle. The reader might observe that these general comments applied equally to the bulk 
side of the correspondence as well as to the boundary side, and might therefore wonder what 
have we gained by translating one hard problem into another hard problem. The purpose 
of this Review is to indicate what in fact we have gained by using a holographic dual, and 
to outline some of the methods that have been used to obtain further understanding of 
dynamical systems. We will structure our presentation of the various approaches according 
to the severity of time-dependence they can handle, i.e. how strongly out-of-equilibrium 
evolutions they apply to. 

A useful starting point is to consider a well-understood global equilibrium situation, and 
try to understand the response of the system under small deviations from this equilibrium. 
If the amplitude of such deviations is suitably small everywhere, the system may be studied 
using a linear response theory, which will be the focus of §3 and §4. We will briefly review 
the basic concepts in linear response theory in §3 focussing on two main aspects: the use 
of retarded correlation function in equilibrium to extract dynamics in the presence of small 
fluctuations, and the behaviour of fast probes in an equilibrated medium as modeled by 
Langevin dynamics. In §4 we will see how this physics can be mapped into the gravitational 
arena. The linear response regime is in fact simplest to understand from the AdS/CFT 
correspondence, for the computation of correlation functions is the best understood part of 
the AdS/CFT dictionary. We will also describe the behaviour of probes and their stochastic 
dynamics by drawing connection with semi-classical dynamics in black hole backgrounds. 

While surprisingly powerful, linear response theory requires small deviations from global 
equilibrium, leaving more general dynamics inaccessible. To go beyond linear response, we 
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Fig. 1: Rough indication of the regimes of validity of the linear response theory and the fluid/gravity 
correspondence, in the space of perturbations from global thermodynamic equilibrium, labeled 
by the amplitude of perturbations A and the wave number (or frequency) k = 4nfp/^. relative 
to the microscopic scale. We have indicated the relevant sections of the paper where the 
different regimes are discussed from the holographic perspective. 

will first focus on long-wavelength IR physics, where the coarse-grained description of an 
interacting field theory is provided by fluid dynamics. Fluid dynamics can of course describe 
fluids well out of global equilibrium, as long as the fluid variables, such as local temperature, 
fluid velocity, etc., make sense. In particular, the spatial and temporal variations of these 
variables must occur on much longer scales than the microscopic ones, but the amplitudes 
need not be small. Borrowing notation from kinetic theory, let us denote this microscopic 
scale ^mfp; and the typical scale on which the fluid variables vary L. Then the regime in which 
the fluid description is meaningful, or equivalently the regime where the system attains local 
thermal equilibrium everywhere, is given by the condition L ^ ^mfp; we will refer to this 
as the 'long-wavelength' regime. How does the fact that we have a fluid description of a 
given configuration help? As we will see in §5, it conveniently constitutes a large truncation 
of the relevant degrees of freedom describing the system. Using the recently-formulated 
fluid/ gravity correspondence [22], there is a one-to-one mapping between any such solution 
to fluid dynamics and a bulk solution to Einstein's equations describing a large non-uniform 
and dynamically evolving black hole. Note that material covered in sections §4.1-§4.2 and 
§5 has previously been reviewed in [23] and [24, 25], respectively. 

Fig. 1 illustrates the two regimes of validity discussed so far. Linear response theory is 
valid for small deformations from global equilibrium, while fluid/gravity correspondence can 
be used in the long-wavelength regime. The main points to note are that these two regimes 
are a-priori distinct, and that they still leave the more interesting region of large deviations 
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from global equilibrium out of reach. 

The general story of understanding physics out of equilibrium then involves consideration 
of deformations that take us to the non-linear regime in amplitude and frequency. This 
translates to the full-blown study of gravitational dynamics, including the fascinating physics 
of black hole formation. In §6 we will review recent progress in understanding various aspects 
of such analyses. Finally in §7 we will take the opportunity to address interesting questions 
about gravity: what does our understanding of physics of strongly coupled field theories 
teach us about gravitational dynamics? Our focus in this brief section will be to describe 
various attempts in the past decade to extract features of the bulk geometry in terms of field 
theory observables. We conclude with a discussion in §8. 

2 The AdS/CFT dictionary 

Before we delve into the details of physics out of equilibrium, we pause to recall some salient 
facts about the AdS/CFT correspondence. As already mentioned, this is well explained in 
the classic reviews [5, 6], so accordingly we will be brief and simply collect facts that are 
necessary for the current discussion. 

The essence of the AdS/CFT correspondence is that strongly coupled field theory dy- 
namics is recorded in terms of string theory (or classical gravity if the field theory admits 
an appropriate planar limit) with appropriate asymptotically Anti-de Sitter boundary con- 
ditions. Spacetimes which are asymptotically AdS can be thought of as deformations of 
pure AdS by normalizable modes in the supergravity description. The framework is also 
sufficiently rich to allow for deformations of the field theory, for instance by turning on rel- 
evant operators.^ The basic dictionary relates the field theory observables, which are gauge 
invariant operators (local or non-local), to their counter-parts in the string (or gravitational) 
description. For instance, local gauge invariant single-trace operators formed out of the fun- 
damental fields of the gauge theory such as Tr {F'^'' F^^) map to single-particle states in the 
bulk spacetime, while non-local operators such as Wilson loops map to string or D-brane 
world-sheets. 

In this language, pure AdS spacetime characterizes the UV fixed point of a quantum field 
theory. The fixed points of interest are field theories with at least M = 1 superconformal 
symmetry.^ The central charge of this CFT is given in terms of the geometry of the AdS 

^ This type of deformation can be used to study field theories on curved backgrounds as described in 
detail in [26] (see also [27, 28, 29] for recent investigations). However, from the bulk perspective, these are 
large deformations, involving asymptotically "locally AdS" geometries. 

^ While non-supersymmctric AdS compactifications also arise in string theory, they are seldom sta- 
ble. Some notable exceptions which satisfy perturbative stability conditions, i.e. spectrum statisfying the 
Breitenlohner-Freedman bound [30], were described in [31]; however to the best of our knowledge non- 
perturbative stability of such vacua have not been established (for instance note that non-supersymmetric 
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spacetime. In the familiar examples arising from string theory one typically encounters 
spacetimes of the form AdSd+ixA" supported by various fluxes, where X is generically a 
compact manifold (which is (9 — (i)-dimensional in string theory or analogously (10 — d)- 
dimensional for M-theory solutions). 

The field theories are characterized by a dimensionless coupling constant (s) and the gauge 
group. We will use the 't Hooft coupling parameter A, while the information regarding the 
rank of the gauge group is given by the central charge c of the CFT. The latter is in turn 
given in terms of the volume of X. Schematically, 

^ / .... ^ . ■ ■ ■ ) ^ ^ / v-i . (2.1) 

Denoting the AdS^+i curvature scale by -RAds and taking into account the basic relation 
between the 10-dimensional Planck and string length scales [33] 

16nG^^'^ = {2nyg'Jl = i2nyil, (2.2) 

we obtain 

^Ads _ Vol(A') Rj^^g 



oc c . (2.3) 



For instance, in the celebrated duality between SU{N) M = 4 SYM and Type IIB string 
theory on AdSs x one finds 

^ -^AdS _ -^AdS _ ^ /9 /|^ 

""^^ 16vrG(5)"(2vr)^^?.^^r Svr^- ^''"^ 

As described above, deformations away from the UV fixed point correspond to deforma- 
tions from pure AdS spacetime. Normalizable modes in the bulk spacetime would be related 
to deformations that are engineered by giving vacuum expectation values to the dual field 
theory operator. Non-normalizable modes in the bulk are non-fluctuating sources that can 
be used to deform the fleld theory Lagrangian.^ Finally, in our considerations we will not 
allow irrelevant deformations of the CFT, as these would correspond to destroying the AdS 
asymptotics. 



2.1 Regimes of interest 



The AdS/CFT correspondence is a profound correspondence between two quantum theories; 
but for general values of the parameters these theories are comphcated and beyond compu- 
tational control. The set of limits we focus on to gain control and insight is the following. 

quotients of AdSsxS^/Zfe do suffer from non-perturbative instabilities [32]). 

^ In certain situations one encounters a choice of boundary conditions to impose on bulk fields. A case 
in point is that of massive scalar fields with mass lying close to the Breitenlohner-Freedman bound [30, 34]. 
We assume that a particular choice has been made among the various possibilities and refer to normalizable 
and non-normalizable with respect to this choice of boundary conditions. 
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• c — )■ oo or equivalently — oo: quantum corrections are suppressed and the gravi- 
tational theory becomes classical. This also has the added merit of suppressing string 
interactions since (7^ ~ 1/c in the planar limit. 

• A — > 00: stringy {a') corrections are suppressed, so the bulk theory is simply Einstein 
gravity interacting with other fields, while the boundary theory describes dynamics 
of local single trace operators. At this level, the field theory dynamics still looks 
complicated and highly nonlocal: the single trace operator expectation value is related 
to the asymptotic fall-off of the corresponding bulk field. Thus although the bulk fields 
evolve according to the 2-derivative field equations, the dynamics on the boundary 
contains infinite number of (spatial and temporal) derivatives. 

Hence in the large- A, large-c regime we are essentially down to studying classical gravi- 
tational dynamics in order to elucidate aspects of the field theory at strong coupling. This 
is of course a drastic simplification of the problem, but nevertheless it warrants further sim- 
plification to gain tractability. Classical gravitational dynamics in asymptotically AdS^i+i 
spacetime involves an infinite number of fields arising from Kaluza-Klein (KK) modes on 
the compact space X. Fortunately, further truncation can often be achieved by the magic 
of consistent truncation [35]. The basic philosophy behind consistent truncations is to find 
an appropriate ansatz for Type II or M-theory fields, which can acturally be reproduced by 
examining the dynamics of a truncated set of fields in d + 1 dimensions. The most familiar 
example of consistent truncation are the gauged supergravity theories which keep only the 
lightest modes under the KK reduction. More complicated examples including massive fields 
have been constructed recently [36, 37] and play an important role in the applied AdS/CFT 
correspondence. 

One hassle with consistent truncations is that the lower-dimensional theory, and hence 
the dynamics of the dual field theory, are model dependent, i.e., dependent on the choice of 
the internal manifold X (and in general on the fiuxes turned on). Although this prevents 
one from making universal statements, it has the opposite advantage of incorporating richer 
dynamics. As we describe later, such extensions are in fact necessary in order to study the 
behaviour of field theories in grand canonical ensembles with prescribed chemical potentials. 
Nevertheless, it is useful to ask whether one can further distill the essential features of the 
correspondence to a minimal classical gravitational Lagrangian, and use this to identify 
universal features that are shared by a wide variety of QFTs. This is in fact possible and 
is achieved in the simplest imaginable manner: by studying the dynamics of pure gravity in 
AdS spacetimes. Let us pause to review that argument before proceeding. 

Given any solution of the form AdS^+i x X one can in fact argue that there is an universal 
sub-sector which simply comprises pure gravitational dynamics in AdS^+i, as can be seen by 
Kaluza Klein reducing on the compact space X. Clearly any non-trivial solution obtained in 
d+1 dimensions uplifts to a solution of the full string theory equations of motion. As a result. 
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we always have a consistent truncation of the comphcated lower dimensional Lagrangian to 
just Einstein gravity with negative cosmo logical constant, where the only dynamical mode 
is the graviton. Since in the field theory dual this graviton mode corresponds to the stress 
tensor expectation value, this truncation provides a decoupled sector with universal dynamics 
for the stress tensor. In other words, the stress tensor obeys the same equations of motion in 
each of such infinite class of strongly coupled field theories. From the field theory perspective, 
the different theories are simply characterized by their differing central charge. Hence apart 
from an overall normalization we will usually be probing dynamics across a wide class of 
theories. 

The complete dynamics of the stress tensor, which in particular allows one to compute all 
the n-point functions of T^j^ in the prescribed state, is still a complicated and non-local system 
from the field theory perspective. After all, this dynamics requires one to be able to solve 
for the dynamics of the non-linear Einstein-Hilbert Lagrangian (with negative cosmological 
constant) in d + 1 dimensions. Therefore to simplify things further, in §5 we will take one 
further limit: we will focus on configurations wherein the stress tensor varies sufficiently 
slowly compared to the local equilibration length scale (the long- wavelength limit). Such 
configurations will then be locally thermalized, allowing for a much simpler description in 
terms of fiuid dynamics. In other words, the equations for the stress tensor in this limit 
reduce to generalized Navier-Stokes equations [22]. 

Once one has an understanding of the dynamics of the field theory stress tensor, or 
equivalently the bulk dynamics of pure gravitational degrees of freedom, one can enlarge the 
system to include other fields. For instance, a natural extension of the canonical ensemble 
in field theories with conserved charges is to include chemical potentials for the said charges 
and examine the behaviour of the grand canonical ensemble. Since global symmetries in 
the field theory map to gauge fields in the bulk spacetime, one naturally ends up studying 
the behaviour of Einstein-Maxell or Einstein- Yang-Mills type theories (again with negative 
cosmological constant of course). In this context one can furthermore consider the behaviour 
of composite gauge invariant operators of the field theory carrying various quantum numbers 
in these ensembles. The gravitational problem then generalizes appropriately to the physics 
of the dual bulk fields. This general scheme has recently been applied to study the phase 
structure of field theories at finite temperature and density (i.e., non-zero chemical potential) 
and has unearthed many interesting features which share qualitative similarities with the 
dynamics of superconducting phase transitions, non- Fermi liquid behaviour etc., which has 
been reviewed in [12, 13, 7, 14]. 

Our strategy in the following will be to consider the simplest setting of just gravitational 
physics in the bulk, since as argued above this sector is universal across a wide variety of 
field theories. To this end, our bulk analysis will involve the dynamics of Einstein gravity 
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with a negative cosmological constant, i.e., 



i^bulk — - 

IGvr Gj^ 

Einstein's equations are given by^ 



I d'+'xV^{R-2A) . (2.5) 



1^^ d{d-^ 

(2.6) 



Emn = Rmn ~ 2^ Gmn ^-j^ Gmn — 



d 

^AdS 



Rmn + — Gmn — . 



If the bulk AdS spacetime geometry is some negatively curved {d + l)-dimensional 
Lorentzian manifold, A4d+i, with conformal boundary dA4d, then the field theory lives on a 
spacetime Bd of dimension d in the same conformal class as dM.d- Choosing an appropriate 
conformal frame, one may identify Bd and dM.d and speak of the field theory as living on the 
AdS boundary. From the standpoint of the bulk theory, the choice of metric on Bd fixes a 
boundary condition that the bulk solution must satisfy.^ The correspondence is simplest to 
state for conformal field theories in dimensions where the trace anomaly vanishes, but with 
appropriate care, the correspondence also holds in the presence of a trace anomaly, and it 
can accommodate non-conformal deformations. 

In general, one could have multiple bulk spacetimes M.d+i whose boundary is the space- 
time Bd on which our field theory lives. In such cases, the AdS/CFT prescription of [4, 38] 
requires that one view all such possibilities as saddle points for the string theory (or gravity) 
path integral and one is instructed to sum over all such possibilities. Of course, the saddles 
might exchange dominance as one changes the boundary manifold; this can be viewed as 
a phase transition of the field theory as one dials an external parameter (in this case the 
geometry of the non-dynamical spacetime which it lives on). We will shortly encounter an 
example of such phase transitions for field theories on compact spatial volume. 



2.2 Field theories in the canonical ensemble & thermal equilibrium 

Finite temperature physics in the field theory can be realized by coupling the system to a 
heat bath, or more precisely by looking at the thermal density matrix 

p = e-^^ (2.7) 

* We use upper case Latin indices {M, A^, • • • } to denote bulk directions, while lower case Greek indices 
{/i, V,- ■ ■} will refer to field theory (or "boundary") directions. Finally, we will use lower case Latin indices 
{«, J, • • • } to denote the spatial directions in the boundary. 

^ In standard AdS/CFT parlance, this amounts to fixing the non-normalizable mode of the bulk graviton 
to obtain the desired metric on Bd- 
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where H is the field theory Hamiltonian. The dual spacetime should have a natural thermal 
interpretation. It is a well-known fact going back to the seminal works of Bekenstein and 
Hawking [39, 40] that black hole spacetimes with non-degenerate event horizons naturally 
exhibit features associated with thermal physics; thereby one is led to expect that black hole 
spacetimes to play a role in describing the dual of a finite temperature field theory, which 
is further supported by the intuition mentioned in §1 that endpoints of generic evolutions 
should match in the dual descriptions. It is, however, logically possible that one also has 
to consider 'thermal geometries' (such as thermal AdS) which just have the Euclidean time 
circle periodically identified. 

To understand this issue better it is useful to think of the thermal density matrix by 
working in Euclidean time. On the field theory side one can achieve finite temperature 
by putting the theory on iS^ = x Ad-i where the Euclidean time circle has period 
(3 = 1/T and Ad^i is the spatial manifold. For compact Ad-i, such as for Ad-i = S'^^^, 
one does indeed have two candidate bulk spacetimes satisfying the boundary conditions [38] . 
The first is the so-called thermal AdS spacetime which is AdS^+i with periodically identified 
Euclidean time coordinate. 



ds^ = - (1 + dr + ^ + dni_-^ , te = -it , te^te + 13 . {2. 

The other saddle point is the static, spherically symmetric Schwarzschild-AdS^+i spacetime, 

2 d-2 / 2 \ 

ds' = -f,{r)dt' + — + r^dnl,, hir) = l + -^--±^(l + -^) , (2.9) 
Jg\' ) -"-AdS ' V -"-AdS/ 

for which the Hawking temperature T is related to the size of the horizon r_)_, as 

drl + (d-2)Rl,^ 
T=^ V J AdS _ 2.10 

These two geometries, (2.8) and (2.9), exchange dominance at r+ = -RAds [38]: at low 
temperature, the thermal ensemble is dual to the thermal AdS spacetime and has a free 
energy of 0{1), while at high temperature, the correct dual is the Schwarzschild-AdS(i+i 
which has a free energy of 0{c). This phase transition is referred to as the Hawking-Page 
transition [41] and is best thought of as a confinement-deconfinement transition (since the 
jump in the free energy is large at large central charge c as required for the planar limit). 
There is furthermore strong evidence that the transition persists to weak coupling where it 
has been identified Hagedorn transition [42]. 

For most of our discussion, however, we are going to be interested in the dynamics of 
field theories on non-compact spacetimes; our focus will typically be on Minkowski spacetime 
^Yiis case there is no phase transition: the fiat space limit is essentially the same 
as the high temperature limit for conformal field theories. As a result, the relevant geometry 
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dual to the thermal density matrix (2.7) in the field theory is the planar Schwarzschild- 
AdSd+i black hole (where WLOG^° we have fixed -RacIS = 1 and r+ = 1): 

ds^ = I -fir) de + J2idxrj + , fir) = 1 - ^ (2.11) 

The causal structure of this solution is easily determined: the spacetime has a spacelike 
curvature singularity at r = 0, cloaked by a regular event horizon at r = 1, and a timelike 
boundary at r = oo. This simple solution is of course static and translationally invariant in 
the boundary directions parameterized by (t,x*). One can in fact generate a rf-parameter 
family of solutions by boosting in R"^"^'^ with normalized d-velocity and scaling r. This 
generates stationary black holes whose horizon size is given (after scaling) by r_|_ and the 
boost velocity u'^ enters into the Killing generator of the horizon via 

r = «'(J,J. (2.12) 

It will turn out that a more convenient form of the metric is one which is manifestly 
regular on the horizon as well as being boundary-covariant. We can obtain such a form by 
starting from (2.11), then change to ingoing Eddington coordinates to avoid the coordinate 
singularity on the horizon: v = t + r^: where dr^: = ;:2y^, and finally 'covariantize' by 
boosting: v — )■ u^x^, Xi — )■ Pi^x^, where P^^ is the spatial projector, P^^ = r]^^ + u^Uu. 
This leads to the form of the metric we will use later: 

ds^ = -2u^,dx^dr + r'^ [r/^^ + (1 - /(r/r+)) u^u^] dx>' dx" , (2.13) 

The event horizon is now at r = r+, which in turn is related to the temperature via the large 
r+ limit of (2.10), 

T=^r,. (2.14) 

Once the bulk black hole solution is determined, it is straightforward to use the holo- 
graphic prescription of [43, 44] to compute the boundary stress tensor. To perform the 
computation we regulate the asymptotically AdS^+i spacetime at some cut-off hypersurface 
r = Ac and consider the induced metric on this surface, which (up to a scale factor involving 
Ac) is our boundary metric Qf^^. The holographic stress tensor is given in terms of the extrin- 
sic curvature K^^, and metric data of this cut-off hypersurface. Denoting the unit outward 
normal to the surface by we have 

K^u = g^pV''n, (2.15) 

Note that in addition to the choice of AdS length scale, the presence of full SO{d — 1, 1) invariance of 
the solution combined with the underlying SO{d, 2) isometry of AdS^+i allows us to rescale the coordinates 
and set the horizon to be located at r = 1. 
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For example, for asymptotically AdS^+i spacetimes, the prescription of [44] gives 



A 



d-2 



lim 



(d+1) 
N 



K"" - K gf"" - {d - 1) g 



1 



d-2 



R^"" --Rg"'' 
2 ^ 



(2.16) 



where K'^'^ is the extrinsic curvature of the boundary. Implementing this procedure for the 
metric (2.13) we learn that the AdS/CFT correspondence maps this bulk solution to an ideal 
fluid characterized by temperature T and fluid velocity u^^. In particular, the induced stress 
tensor on the boundary is 



(2.17) 



where Ad is an overall normalization that is proportional to the central charge c of the fleld 
theory. Note that this stress tensor is traceless, T'^ = 0, as indeed is expected for a CFT. 



2.3 Equilibrium in grand canonical ensembles 



A natural extension of the above framework is to consider systems at finite density by 
generalizing (2.7) to 

p = ^-P(H-^rQ') (2.18) 

where fii are chemical potentials for the conserved charges . The choice of chemical 
potentials is of course restricted by the global symmetries of the field theory. If we insist 
on restricting attention to the universal dynamics of the stress tensor alone, then the only 
possible generalization is to consider chemical potentials for rotations Qj. For a field theory 
on Ad-i = S'^~^ one in general has independent rotations and each of these can be 

given a non-zero angular velocity Qj. 

However, if we are willing to add matter fields in the bulk, i.e. extend (2.5) to include 
additional (non-gravitational) degrees of freedom, then we can generalize the discussion to 
include more interesting chemical potentials. First of all, we note that in order to incorporate 
chemical potentials, the field theory must admit conserved currents. Conserved currents in 
the field theory map to gauge symmetries in the bulk spacetime.^^ For every conserved field 
theory current Jj one therefore has a bulk gauge field Aj^. The bulk gauge field obeys 
some equations of motion in the AdS^+i spacetime and its non-normalizable mode in the 
near-boundary expansion corresponds to the boundary chemical potential fij. 

Typically, large-A^ field theories with holographic duals have conserved charges; for in- 
stance M = A SYM has a 5*0(6) ^ SU (4) global symmetry and one can consider the grand 
canonical ensemble with /^/(s) corresponding to Cartan subgroup of this non-abelian global 

^-'^ This is in fact reminiscent of the theorem "No global symmetries in string theory", which follows from 
the observation that conserved currents in spacetime lead to world-sheet vertex operators of weight (1,0) 
and (0, 1) respectively, which in turn imply massless particle states in spacetime. 
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symmetry group. In fact, for a wide variety of A/" = 1 superconformal field theories in d = 4 
one has three conserved U{1) charges which geometrically can be related to the fact that 
these field theories arise from compactifications of Type IIB supergravity on toric Sasaki- 
Einstein manifolds. A simple example to keep in mind is the Einstein-Maxwell Lagrangian 
(with perhaps Chern-Simons terms in odd bulk dimensions) which can be used to study the 
grand canonical ensemble with U{1) charge chemical potentials, which can be obtained via 
a consistent truncation of gauged supergravity theories. 

Given an appropriate definition of the field theory grand canonical ensemble, the equi- 
librium solution corresponding to this ensemble can be found by looking for static black 
holes carrying the appropriate charges. For Einstein-Maxwell theory these would be just the 
Reissner-Nordstrom-AdS (RNAdS) black holes, and the thermodynamic properties of these 
black holes correspond to the equilibrium thermodynamics of the field theory. 

In the presence of matter, one could have interesting phase transitions over and above the 
analog of the Hawking-Page transition discussed in §2.2. For pure Maxwell field interacting 
with gravity, these include charge redistribution as discussed originally in [45, 46], or the more 
recent examples involving Chern-Simons dynamics [47]. Inclusion of charged scalars (and 
sometimes neutral scalars) can also lead to interesting phase transitions as originally pointed 
out by [48], which plays an important role in the physics of holographic superconductors 
[49, 50]. 

3 Linear response theory 

Having reviewed the basic framework of the AdS/CFT correspondence we now turn to de- 
scribing the essential features of linear response theory. This will set the stage for our 
discussion of deviations away from equilibrium in §4. 

Consider a generic quantum system characterized by a unitary Hamiltonian H acting on 
a Hilbert space "H in equilibrium. The system could be in a pure state | G "H or more 
generally in a density matrix p; the only requirement is that the system be in a stationary 
state. We now wish to perturb the system away from equilibrium and analyze the dynamics. 
A general deformation can be thought of as a change in the evolution operator H — )■ H+Hpert 
and one is left with having to examine the dynamics with respect to this new Hamiltonian. 
Things however simplify if we can focus on deviations which are small in amplitude - this is 
the regime of linear response theory, which we now briefly review. For a beautiful account 
see the original derivation by Kubo [51, 52] and the review [53]. 

In the linear response regime, one imagines the system being perturbed by a weak external 
force. To wit, one can write -ffpert = —AJ^{t) where J^{t) is the external force (with time 
dependence explicitly indicated) and A is the canonical conjugate operator to the force. In 
the linear response regime the amplitude of the force is constrained to be small, |v^(t)| — 
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^ ^ 1, while we allow arbitrary temporal dependence, which in particular could involve 
high frequency modes being excited; see Fig. 1. We essentially wish to do time dependent 
perturbation theory (being careful of causality in relativistic theories) for such deviations 
away from equilibrium. 



3.1 Response functions for deviations from equilibrium 

To monitor the departure from the stationary state, we can pick some observable, say the 
expectation value of a local operator O in the theory. For concreteness, we will assume 
that the system was in a density matrix p before we perturbed it. Stationarity of this state 
demands that [H, p] = 0. Turning on the perturbation //port will cause a deformation to the 
density matrix. Suppose that we encounter a new density matrix p, which now has to satisfy 
Hamiltonian evolution with respect to H + //pert, i-e., 

Z^p=[// + //pert,p] (3.1) 

For small amplitude perturbations, we can assume that p = p + Sp and solve for 6p formally 
in terms of a time ordered integral expression: 

5p{t)=i f (it'e-^(*-*')^ [Ap] e'(*-*')^J^(t') (3.2) 

J —oo 

which is derived by solving the linearized version of the evolution equation (3.1). 

Once we have the change in the density matrix, it is easy to compute the change in any 
measurement occurring due to the perturbation. We can estimate the response of the system 
by examining the expectation value of some observable O, which can be obtained directly 
as: 

(0(t)) = Tr(pO) 

^(5(C(t)) =iTr dt'[A,p]0{t-t')J^{t')^ (3.3) 

It is conventional to define the response function TZ (sometimes called the after-effect func- 
tion) as the change in the operator expectation value for a delta function perturbation; i.e., 
for F{t) = S{t). From (3.3) one recovers 

7^^(t) = t Tr {[A, p] 0{t)) ^ -I Tr (p [A, 0{t)]) (3.4) 

where we have tried to make clear in the notation the idea that one measures the response 
of O to a perturbation caused by the deformation due to a force conjugate to A. By taking 

One could also consider non-local operators and more exotic observables such as entanglement entropy; 
we will discuss the latter in §7. 
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a Fourier transform of the response function, one arrives at the admittance: 

POO 

x^M = hm / 7^^(^)e-*^*-^*c^^ (3.5) 

It should be clear from the above discussion that this formalism is sufficiently general to 
accommodate arbitrary linear changes in a given dynamical system. Non-linear corrections 
can be explored in perturbation theory, generalizing familiar ideas from time dependent 
perturbation theory in quantum mechanics. 

3.2 Retarded correlators & Kubo formulae 

The perturbations discussed so far are explicit perturbations on the system caused by some 
external force J-'. One can as well envisage a perturbation driven purely by thermal fluc- 
tuations, which are not a-priori related in any obvious way to external forces acting on 
the system. However, thermal fluctuations can be measured by looking at the response of 
the extensive thermodynamic variables to variations in the local energy or charge densities. 
These in fact can also be treated in linear response and lead to the famous Kubo formulae 
for the transport coefficients. We now give a brief overview of these concepts valid in any 
field theory; in §4 we will demonstrate how these formulae can be applied in the AdS/CFT 
context. 

Let us return to the response function given in (3.4) and extract some essential features 
that we will use in later analysis. While that discussion was sufficiently general, it is worth 
translating this into more familiar language. Physically we wish to monitor the behaviour of 
the expectation value of some operator Oa when the system is subject to some perturbation. 
We can imagine this occurring via a direct coupling of the operator to some sources J'a, 
whose effect is to change the action via: 

S = So + j d^xJa{x)Oa{x) (3.6) 

The response of the system due to this change can be obtained by rewriting the result (3.4) 
slightly. Since the system will respond only after the perturbation, causality demands that 
we simply convolve the retarded Green's function of the operator Oa to the sources Jb that 
causes it to deviate from stationarity. In particular, we have 

{Oa{x)) = - j d'yGl{x,y)My) (3.7) 

where the retarded correlation function is defined as usual via 

Gf,(x, y) = -I e{x' - y') { [a(x), O^iy)] ) (3.8) 

Therefore given the retarded correlation functions, one can immediately infer, in the linear 
response regime, the manner in which the system under consideration reacts to the external 
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disturbance caused by the sources. One will often be interested in the behaviour of the 
Green's functions in momentum space. To this end we define the Fourier transform of the 
retarded correlators via 



where we have assumed translational invariance and defined the d- vector fc^ = {u, k). 

A class of observables that are interesting to examine are those corresponding to conserved 
currents in the theory, viz., the energy-momentum tensor T^^, or generic conserved global 
charges J^. If we consider systems in thermal equilibrium, where deviations from thermality 
are engineered by density or charge fluctuations, then one is naturally led to studying the 
retarded Green's functions of these conserved currents: 



G^u,ap{x, y) = -i Tr (p [T^^,{x),To,p{y)]) = -i9{x -y){ [r^^(x), Tc,^(?/)] )/? 
Gf,y^a{,x,y) = -iTi{p[T^^{x),Jaiy)]) = -i9ix^ -y^){[T^uix),Jaiy)])(3 

G^,aix,y) = -tTTip[Ux),Uy)]) = -i0ix'-y'){[U^),Uy)])f3 (s.io) 



where we have assumed that the density matrix in question is appropriate to the ensemble 
under consideration. The correlators indicated on the RHS in (3.10) are therefore the thermal 
correlators, with perhaps fixed chemical potentials. 

Interesting physical quantities characterizing the system can be extracted by examining 
the momentum dependence of the retarded correlation functions of the conserved currents. 
The momentum space correlators have non-trivial analytic behaviour, with poles in the 
complex k plane. One can read off the dispersion relation for the associated modes, by 
solving for Ci;(k). Since thermodynamic systems typically incorporate dissipative effects, 
these dispersion relations typically have imaginary pieces which capture the rate at which 
the system relaxes back to equilibrium. Stability of the quantum system demands that the 
perturbations damp out exponentially in time. This translates to the poles of the retarded 
Green's functions lying entirely in the lower half plane of the complex frequency space. We 
will see shortly that these poles of the retarded Green's functions are in fact associated with 
the quasinormal modes of black hole geometries (which describe how a black hole settles 
back to its quiescent equilibrium state) via the AdS/CFT correspondence. 

To be specific, let us consider a four dimensional conformal field theory and examine the 
stress tensor retarded correlator. It is useful to take into account the stress tensor conserva- 
tion equation which implies a Ward identity k'^^G^p^ais = 0.^'^ It is in fact convenient to pick 
a direction in momentum space, say k = ke^, and describe modes as being longitudinal or 
transverse to this choice of momentum. One can then show that the transverse components 

We are here ignoring contact terms which can appear in the Ward identities. In momentum space these 
wiU give rise to analytic pieces and will therefore be irrelevant to our discussion of of poles. 




(3.9) 
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of the stress tensors have the behaviour 



{u, k)} 



- j/;;^, —u k, uj^ 
\^Qz{^,k) 



} 



Gi{(^,k) 
ij^ — k'^ 




(3.11) 



while the longitudinal modes behave as 



{Gu,tt{oJ, k), Gtt,tz{(^, k), Gtz,tz{uj, k)} 



- {A;^, —CO k, up' 



} 



k^ 02(00, k) 
(a;2 - k^y 



(3.12) 



The correlators are thus completely described by three scalar functions of momenta, Qi{k) 
for z = 1, 2, 3, which exhibit the aforementioned poles. 

Out of the set of poles of the retarded Green's function, of special interest are those 
that capture the late time behaviour. These necessarily involve small imaginary parts in 
Ci;(k) (since the modes with large imaginary parts damp out quickly and therefore have short 
half-lives). These special set of poles are the hydrodynamic poles; they capture, and in fact 
provide a basis for, a complete description of the interacting quantum system via linearized 
hydrodynamics.^^ In the low-frequency regime, i.e., for w, <^ T, it turns out that the 
function Q^lujk) defined in (3.11) is non-singular, while the other correlation functions 
exhibit poles. In fact, at low frequencies only the pieces of the correlation functions that 
correspond to conserved quantities can be singular, for it is these modes which due to the 
conservation law take a long time to relax back to equilibrium. 

The hydrodynamic poles of the system are characterized by having dispersion relations 
which satisfy the constraint uj{k) ^ as k ^ 0. Fluctuations transverse to the direction of 
momentum flow, captured by Qi{uj, k), give rise to dispersion relation of the form 



which is characteristic of a diffusive mode, with diffusion constant D. For energy- momentum 
transport this is the shear mode, with the diffusion constant D being related to the shear 
viscosity rj of the system via D = , where e and P are the equilibrium energy density and 
pressure, respectively. The longitudinal component, given by ^2(1^, k), has poles at locations 



which describe sound propagation in the medium with velocity Vs and attenuation F^. Note 
that the sound mode is the only propagating mode in the hydrodynamic limit. 

While one can examine the analytic structure of the retarded correlators in momentum 
space and extract interesting transport properties, it is useful to obtain direct formulae for 

Note that hydrodynamics is a good approximation to any physical system close to equilibrium, for 
fluctuations that are sufficiently long in wavelength. We discuss this in detail in §5. 



Diffusive mode: u = —i D k'^, 



(3.13) 



Sound mode: 



u = ±f , k — iV^k"^ 



(3.14) 
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them. These are the famous Kubo formulae. For instance to find the shear viscosity of the 
system, one simply takes an appropriate zero frequency limit of the retarded correlator, i.e., 

r/ = -lim -^lm(G5^,,^(a;,0)). (3.15) 
Similar expressions can be written down for the charge conductivity, etc.. 



3.3 Brownian motion of probes and Langevin dynamics 

Our discussion thus far has been anti-chronological in a historical sense, for we have used 
general notions of quantum field theories and statistical mechanics to arrive at the response 
of the system to external perturbations. Historically, these ideas were first explored in kinetic 
theory, where it was realized that one could systematically account for the deviations away 
from purely thermal behaviour. We will refrain from repeating these ideas here in the context 
of departures from equilibrium of the system as a whole, but instead use these concepts to 
describe the physics of a single probe particle in a thermal medium. We have in mind a 
projectile that moves through some plasma medium. The medium itself will be taken to be 
in a thermal ensemble and we will be interested in the manner in which the probe particle 
loses energy to the medium. Moreover, it is also well known that even if the particle attains 
equilibrium with the plasma, it will continue to be buffeted by thermal fluctuations from 
the medium and undergo random motion. This is the famous stochastic Brownian motion, 
which is best described in the limit of a heavy probe particle in a thermal medium. 

Let us therefore consider the Langevin equation, which is the simplest model describing 
a non-relativistic Brownian particle of mass m in one spatial dimension: 

pit) = -^oPit) + Rit) . (3.16) 

Here p = mx is the (non-relativistic) momentum of the Brownian particle at position x and 
time t, and ' = d/dt. The two terms on the right hand side of (3.16) represent friction and 
random force, respectively, and 70 is a constant called the friction coefficient. One can think 
of the particle as losing energy to the medium due to friction, and concurrently getting a 
random kick from the thermal bath, modeled by the random force. We assume the latter to 
be simply white noise with 

(i?(t)) = 0, {R{t)R{t')) = KoS{t~t'), (3.17) 

where Kq is a constant. Note that the separation of the force into frictional and random parts 
on the right hand side of (3.16) is merely a phenomenological simplification - microscopically, 
the two forces have the same origin, namely collision with the fluid constituents. 

Assuming equipartition of energy at temperature T = ( m i;^ ) , one can derive the follow- 
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ing time evolution for the square of the displacement s{t) [54]: 



: s{tY ) ^ { [x{t) - x(0)]2 ) = ^ 1 - 1 + 

7o 
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(3.18) 



where the diffusion constant D is related to the friction coefficient 70 by the Sutherland- 
Einstein relation: 

D = — . (3.19) 
7om 

We can see that in the ballistic regime, t <ti I/70, the particle moves inertially {s ~ t) with 
the velocity determined by equipartition, x ~ ^JT jm. On the other hand, in the diffusive 
regime, t 3> I/70, the particle undergoes a random walk (s ~ -\/t). The transition is not 
instantaneous because the Brownian particle must collide with a certain number of fluid 
particles to get substantially diverted from the direction of its initial velocity. The crossover 
time between the two regimes is the relaxation time 

1 . , 

^relax ~ , (3.20) 

7o 

which characterizes the time scale for the Brownian particle to forget its initial velocity and 
thermalize. One can also derive the important relation between the friction coefficient 70 
and the size of the random force kq. 

70 = ^ , (3.21) 
2 m i 

which is the simplest example of the fluctuation-dissipation theorem and arises precisely 
because the frictional and random forces have the same origin. 

In n spatial dimensions, the momentum p and force R in (3.16) are generalized to n- 
component vectors, and (3.17) is then naturally generalized to 

(i?,(t) ) = , {R,{t)Rj{t') ) = Ko5ij5it - t') , (3.22) 

where i, j = 1, . . . ,n. In the diffusive regime, the displacement squared scales as 

{s{tf) ^2nDt . (3.23) 

On the other hand, the Sutherland-Einstein relation (3.19) and the fluctuation-dissipation 
relation (3.21) are independent of n. 

Let us now return to the case with one spatial dimension [n = 1). The Langevin equation 
(3.16), (3.17) captures certain essential features of physics, but nevertheless is too simple to 
describe realistic systems, since it assumes that the friction is instantaneous and that there 
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is no correlation between random forces at different times (3.17). If tlie Brownian particle is 
not infinitely more massive than the fluid particles, these assumptions are no longer valid; 
friction will depend on the past history of the particle, and random forces at different times 
will not be fully independent. We can incorporate these effects by generalizing the simplest 
Langevin equation (3.16) to the so-called generalized Langevin equation [53, 55], 

p{t) = - I dt' 7(t - t') p{t') + R{t) + K{t) . (3.24) 

J — oo 

The friction term now depends on the past trajectory via the memory kernel 7(t), and the 
random force is taken to satisfy 

(i?(t)) = 0, {R{t)R{t')) = K{t-t') , (3.25) 

where K{t) is some function. We have in addition introduced an external force K{t) that can 
be applied to the system. 

Let us briefly indicate how, faced with such a system with only K under our control, we 
would extract the physical information, namely 7(t) and ^(t), as well as the characteristic 
timescales involved in collision and relaxation. The latter will offer more direct insight 
into the nature of the medium under consideration. To analyze the physical content of the 
generalized Langevin equation, it is convenient to first Fourier transform (3.24), obtaining 

. m±J^ . (3.26) 

7[u;J — too 

where p{u), R{u), K{u) are Fourier transforms, e.g., 

/oo 
dtp{t)e'^\ (3.27) 
-oo 

while 7[ci;] is the Fourier-Laplace transform: 

POO 

7[w] = / dt-f{t)e''^' . (3.28) 
Jo 

If we take the statistical average of (3.26), the random force vanishes because of the first 
equation in (3.25), and we obtain 

{p{u) ) = fi{u)K{u), ^ / . , (3.29) 

where /x(w) is called the admittance. The strategy, then, is to first determine the admittance 
/x(a;), and thereby 7[w], by measuring the response (^(w) ) to an external force. For example, 
if the external force is 

K{t) = Ko e"^"*, (3.30) 



24 



then ) is simply 

{p{t))=fi{co)Koe~'^' . (3.31) 
For a quantity 0(t), we define the power spectrum Io{^) by 

/oo 
dt ( O(to) O(to + t) ) e*"* • (3.32) 
'OO 

Note that ( O{to) 0(tQ+t) ) is independent of to in a stationary system. The knowledge of the 
power spectrum is equivalent to that of 2-point function, because of the Wiener-Khintchine 
theorem: 

( 0{u) 0{uj') ) = 27:6(00 + u') Io{uj) . (3.33) 
Now consider the case without an external force, i.e., K = 0. In this case, from (3.26), 

piu) = . (3.34) 
7[a;J — too 

Therefore, the power spectrum of p and that for R are related by 

W = I /f""^. 12 • (3.35) 

Hence, combining (3.31) and (3.35), one can determine both 7(t) and /t(t) appearing in 
the Langevin equation (3.24) and (3.25) separately. However, these two quantities are not 
independent but are related to each other by the fluctuation-dissipation theorem, generalizing 
the relation (3.21), cf., [53]. 

For the generalized Langevin equation, the analog of the relaxation time (3.20) is given 

by 



t 



relax 



dt-f{t) 



-1 



/i(w = 0) . (3.36) 



J l[^ = 0] 

If 7(t) is sharply peaked around t = 0, we can ignore the retarded effect of the friction term 
in (3.24) and write 

/ dt' 7(t - t') p{t') ^ / dt' 7(t') ■ p{t) = p{t) . (3.37) 

Jo Jo t-relax 

The generalized Langevin equation (3.24) then reduces to the simple Langevin equation 
(3.16), so that treiax corrcspouds to the thermalization time for the Brownian particle. 

Another physically relevant time scale, the microscopic (or collision duration) time tcoib 
is defined to be the width of the random force correlator function K{t). Specifically, let us 
define 



^coU 







dt ^ . (3.38) 
k{0) 
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If K{t) = k(0) e~*/*<:°", the right hand side of this precisely gives tcoii- This tcoii characterizes 
the time scale over which the random force is correlated, and thus can be thought of as the 
time elapsed in a single process of scattering. In many cases, it is natural to expect that 

^rolax ^ '^coll ; (3.39) 

since, after all, we indicated that it takes a heavy probe many collisions to thermalize. 
Typical examples for which (3.39) holds include settings where the particle is scattered 
occasionally by dilute scatterers as described by kinetic theory, and settings where a heavy 
particle is hit frequently by much smaller particles [53]. However, as we will discuss in §4.3.2, 
for the Brownian motion dual to AdS black holes, the field theories in question are strongly 
coupled CFTs and in fact (3.39) does not necessarily hold. There is also a third natural time 
scale ^mfp 

given by the typical time elapsed between two collisions. In the kinetic theory, 
this mean free path time is typically between the single-collision and relaxation time scales, 
^cou <^ ^mfp ^ ^reiax; but again, this hierarchy is not expected to hold beyond perturbation 
theory. 

The basic message to take away from this discussion is that the linear response regime 
is accessible once one understands the dynamics in equilibrium. The response functions are 
simply given in terms of Green's functions evaluated in the stationary configuration. 

4 Linear response from AdS/CFT: Probes of thermal plasma 

We now proceed to put together the toolkits we presented in the preceding two sections. In 
§3 we have described the basic methods employed in non-equilibrium statistical mechanics 
to understand the physics of systems out of equilibrium, viz., linear response theory. One 
of the fundamental tenets in this approach is that for small-amplitude deviations, one can 
compute relevant observables by computing appropriate correlation functions in the equilib- 
rium ensemble. Since the time-evolution part of the problem has been effectively dealt with 
in this manner, one is left with a much easier task in general. 

However, the computation of equilibrium correlation functions is not as trivial as it 
sounds, especially in circumstances where the underlying quantum system is intrinsically 
strongly coupled. One therefore requires some further insight to deal with such situations. 
Fortunately, for a class of field theories which have holographic duals, the gauge/gravity 
correspondence comes to rescue, as indicated in §2. In fact, one the earliest developed 
technologies within the correspondence was the recipe to compute correlation functions of 
gauge invariant local operators in field theories using their dual gravity picture. In the 
present context this means that the gauge/gravity correspondence provides an efficient way 
to compute the correlation functions relevant for the linear response theory directly in terms 
of classical computations in an asymptotically AdS spacetime. 
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We will begin by a brief review of the techniques employed to compute correlation func- 
tions in the AdS/CFT correspondence, followed by a discussion of lessons learnt by examining 
the linear response regime. Finally, we will discuss how one can monitor the behaviour of 
probe motion (both ballistic and stochastic) in a strongly coupled plasma medium. 

4.1 Computing correlation functions in AdS/CFT 

Let us consider a local gauge-invariant single trace operator 0{x) with conformal dimension 
A in the boundary CFT. To compute correlation functions of this operator, one would deform 
the CFT action by adding a term J d'^x (f)o{x) 0{x) and obtain the generating function of 
the correlators as a functional of the sources (poi^)- The requirement that the deformation 
term be dimensionless implies that (poi^) has scaling dimension d — A. 

In the AdS/CFT correspondence, a given boundary operator maps to a bulk field whose 
spin s is determined by the Lorentz transformation property of the operator O in question. 
The bulk field (j){x, r), with r being the radial coordinate in AdS^+i, has mass m^(A, s). For 
scalar operators one has the relation [4]^^ 



Similarly, for a p-form operator on the boundary, the relation between the conformal dimen- 
sion of the operator A and the mass of the bulk field m is given as: 



Given this map between fields and operators we can go ahead and use the AdS/CFT 
correspondence to compute the generating function of correlation functions. 



where we have schematically indicated the path integral over the quantum fields of the CFT. 

The statement of the AdS / CFT correspondence asserts that this generating functional is 
given by the partition function ^string of the string theory with the fields 0(r, x) prescribed 
to take on the boundary values (poi^x) at the boundary of AdS^+i. In particular, in the 
limit when classical gravity is a good approximation, the string partition function simply 
reduces to the on-shell action of gravity evaluated on the solution to the field equation. The 
on-shell action is usually divergent since we are turning on mode that is non-normalizable 

This is not required to be true for rn^p < m? < m?gp + 1 where ni^gp = is the Brietenlohner- 

Freedman mass, providing the lower bound on the mass of a scalar field in AdS^+i- In the said range, one can 
equally well associate bulk field of mass m with a boundary operator of dimension 
as discussed in [34], which corresponds to an alternative quantization of fields in AdS^+i [30]. 





( A + p) ( A + p - d) = . 



(4.2) 




(4.3) 
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to act as a source. To ensure that we capture the correct physics, we can regulate the AdS 
spacetime at r = and demand that we satisfy the boundary condition there, i.e., demand 
(/)(r, x) — ?► (f)o{x) as r — > e^^. Thus one arrives at the relation derived in [3, 4]: 

(eW M^)oi.) ^^^^ ^ 2,^^.^^ [0(r = e-\x) = 0o(x)] (4.4) 
or in the limit where classical gravity in the bulk is a good approximation 

W[(f)o] = - logZstring - extremum^(^=,-i)=0„ /sugra(0o) (4.5) 

The general scheme we have outlined above works well for computing Euclidean corre- 
lation functions in asymptotically AdS spacetimes, but there are certain subtleties to keep 
in mind while computing retarded correlation functions. A clear prescription was initially 
given in [56] for which a nice supporting argument based on Schwinger-Keldysh contours 
was provided in [57]. More recently, these arguments have been revisited in [58] and a com- 
pact expression for computing two-point functions was provided in [59, 60]. Formal studies 
of these correlation functions from a holographic renormalization scheme and general con- 
tour prescriptions for higher-point functions were discussed in [61, 62]; recently three-point 
functions at finite temperature were computed in [63]. 

Since we will be primarily interested in addressing issues in linear response theory, let 
us record here the prescription derived in [59], relating the retarded Green's function to 
a simple ratio involving the field and its conjugate momentum. In particular, for massive 
sscalar fields in asymptotically AdS spacetimes, 

n(r, k) 



Gnik) = finite! 1™ r^^^"'^) 



Here 11 is the canonical momentum conjugate to the field (which itself is dual to the 
operator O under consideration) under radial evolution in AdS. Furthermore, 0soi('^, k) is 
the on-shell solution to the appropriate wave equation subject to the boundary conditions 
that it be regular in the interior^^ and approaching some chosen boundary value 0o at the 
boundary of the AdS spacetime. The constraints 0o = (i.e. switching off the source) and 
the limit r — )■ oo are necessary to obtain the boundary observable as one anticipates from 
(4.5). We should note that this formula has been written down after taking into account 
the intricacies of the holographic renormalization and hence one is instructed to extract the 
finite part of the bulk calculation. For details on these techniques we refer the reader to 
[64, 65]. 



4.2 Retarded correlators and black hole quasinormal modes 

Given the utility of the AdS/CFT correspondence in computing correlation functions, let us 
now return to the issue of linear response around a given equilibrium configuration. As we 

In the case of spacetimes with horizons, this amounts to demanding that the field be infaUing at the 
horizon, as we explain in §4.2. 
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have seen in §2.2 and §2.3, physics of thermal equihbrium is captured by stationary black 
hole spacetimes in the dual geometric description. Therefore, linear response behaviour in 
the field theory translates directly to the behaviour of linearized fluctuations of bulk fields 
on AdS black hole backgrounds. 

One of the first steps in this direction was taken in [66], who pointed out the connection 
between AdS black hole quasinormal modes and the rate at which disturbances away from 
equilibrium re-equilibrate. Since this connection underpins much of the linear response 
theory we are about to describe, we will pause to recall the basics of the quasinormal mode 
spectrum in black hole spacetimes. 

Physically, quasinormal modes correspond to the late-time "ringing" of the black hole 
geometry. In particular, perturbations of the black hole undergo damped oscillations, whose 
frequencies and damping times are entirely fixed by the geometry and the nature of the 
propagating field, i.e., the modes are determined by the linearized wave operator and are 
independent of the initial perturbation. In fact, it is well understood that black hole space- 
times, owing to the presence of an event horizon into which the fields can dissipate, act as 
open systems; the corresponding spectrum of fluctuating modes is complex. The reason for 
this behaviour is intuitively easy to understand. In classical general relativity, the event 
horizon acts as a one-way membrane; fields fall into the black hole but do not emerge out. 
Mathematically, this translates to an infalling boundary condition on fields at the horizon 
in this black hole background. These same fields are also required to be normalizable near 
the AdS boundary, for one wishes to retain the AdS asymptotics (and therefore in the field 
theory side retain the UV fixed point CFT unperturbed by relevant or irrelevant opera- 
tors). Quasinormal modes for a classical field $ (suppressing Lorentz indices) are defined as 
eigenf unctions of the linearized fluctuation operator which acts on $ in the black hole back- 
ground, satisfying these boundary conditions i.e, ingoing at the horizon and normalizable at 
infinity.^'' 

As initially described in [66] , the quasinormal modes of AdS black holes capture the rate 
at which the field theory, when perturbed away from thermal equilibrium, returns back to 
the quiescent equilibrium state. In this context, one usually concentrates on the lowest set of 
modes, as these dominate the long-time behaviour. Nevertheless, it is possible to give a clear 
interpretation to the entire quasinormal mode spectrum. As was pointed out in [69] in the 
context of 1-1-1 dimensional boundary CFTs and asymptotically AdSs BTZ black holes, the 
entire quasinormal mode spectrum maps to the poles of the retarded Green's functions of 
operators in the canonical ensemble. This was extended to higher dimensions in the seminal 
work of [56] and further elaborated upon in [70]. 

While the relation between quasinormal modes and poles of retarded Green's functions 

An excellent summary of known spectra of quasinormal modes for various fields in diverse black hole 
backgrounds, together with their implication for both astrophysical black holes and in the AdS/CFT context 
can be found in [67]. See also [68] for an earlier discussion of black hole quasinormal modes. 
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holds in general for any operator in the dual field theory, it takes on interesting hues for 
the case when the dual operator corresponds to a conserved current. As discussed in §3.2, 
the analytic structure in the retarded Green's functions of the stress tensor which corre- 
sponds to the hydrodynamic modes of the system, has complex dispersion relations u}{k) 
characterized by the long-wavelength behaviour uj{k) — )■ as A; — )■ 0. This was explicitly 
verified by the computation of gravitational quasinormal modes in planar AdS^+i black hole 
backgrounds, which have translationally invariant horizons and allow for arbitrarily long- 
wavelength modes. On the contrary, global AdS^+i black holes have horizons of spherical 
topology and correspond to field theories living in finite volume on ESU^ = R x S'^^^. One 
then encounters IR effects coming from the finiteness of spatial volume which precludes the 
existence of quasinormal modes with vanishing frequencies. In order to see the hydrody- 
namic behaviour one has to scale the curvature to zero as well, which reduces the problem to 
the planar case; one can systematically account for the curvature corrections as we indicate 
in §5. 

The fact that black hole quasinormal mode spectrum admits modes with hydrodynamic 
dispersion relation leads one to suspect that one can use the gravity analysis to compute 
properties of the fluid description. Indeed as we have sketched in §3.2, one can use the 
behaviour of the retarded Green's functions at zero momentum to learn about transport 
coefficients like viscosity. This was flrst carried out in the AdS/CFT context in [71, 72]. The 
analysis was ground-breaking in that it not only verifled the general intuition that one can 
relate the classical dynamics in a black hole background to the physics of strongly coupled 
plasmas, but it also paved the way for what is perhaps the most famous conjecture in the 
subject, viz., the bound on the ratio of shear viscosity to entropy density, rj/s > [15]. 
For a large class of two-derivative theories of gravity one flnds by direct computation that 
this bound is in fact saturated, which prompted [15]. Understanding its implications and its 
raison d'etre has been the focus of a large body of literature, which we cannot do justice to 
here and point the reader to the excellent review article [23] for developments till a couple 
of years ago. 

This rather small value of shear viscosity obtained in the holographic computations has 
been instrumental in forging connections with ongoing experimental efforts to understand 
the state of matter, the quark-gluon plasma (QGP), produced in heavy-ion collisions and 
RHIC and soon at LHC. Fits to data from the STAR detector at RHIC suggests that the 
QGP behaves close to the deconflnement transition in QCD as a nearly-ideal fluid with very 
low viscosity (see [73] for a discussion of near perfect fluidity in physical systems). This has 
prompted a concentrated effort in the literature and spurred the growth of the AdS/QCD 
enterprise; we refer the reader to the reviews [9, 11] for these developments. 

Recently, this bound has been shown to be violated in higher-derivative theories of grav- 
ity: there are example toy-models such as Gauss-Bonnet gravity [74, 75] and other higher- 
derivative theories [76, 77, 78, 79, 80, 81] and also some string theory inspired constructions 
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of large- superconformal theories [82, 17]. The general consensus at the stage of writing 
this review seems to be that the bound, whilst robust in the two-derivative approximation 
(which corresponds to the strong coupling, large- theory), could in general be violated by 
finite- (string interactions) and also perhaps by a' (large finite coupling) effects. 

The quasinormal mode analysis can also be used to go to higher orders in the hydro- 
dynamic expansion. After all, for a planar black hole one has a spectrum of poles of the 
retarded Green's function which can be used to extract an exact non-linear dispersion re- 
lation u{k) beyond the leading long-wavelength A; — )• approximation. Such techniques 
were first employed in the analysis of [83] who computed certain second-oder transport co- 
efficients. Similar analyses were also carried out in [84, 85]. In an nice calculation [86] 
described the behaviour of low lying quasinormal modes (the modes closest to the real axis) 
and displayed how it exchanges dominance with the next quasinormal mode at some finite 
value of momentum. This in particular indicates the regime of validity of the linear hydro- 
dynamic approximation; for the higher quasinormal modes, while still giving poles of the 
retarded Green's functions, are not part of the effective hydrodynamic theory. 

In summary, there is a direct relation between the physics of black hole quasinormal 
modes and the retarded Green's functions of local gauge invariant operators. For a given 
operator O on the boundary, one identifies the corresponding bulk field (p and computes 
its quasinormal mode spectrum to infer the location of the poles of the retarded Green's 
function. From the retarded Green's functions of conserved currents one learns that the long- 
wavelength behaviour of the interacting CFTs which fall within the purview of the AdS/CFT 
correspondence is described by linearized hydrodynamics. Thus black hole quasinormal 
modes provide a powerful computational technique to learn about the dynamics of strongly 
coupled gauge theories and their relaxation back to equilibrium. They also confirm the 
intuition that the thermal behaviour of strongly coupled field theories can be captured in 
effective field theory by viewing the system as a plasma medium. 

4.3 Probes in the plasma: Dissipation and stochastic motion 

Thus far, we have discussed the behaviour of retarded correlation functions of local, gauge 
invariant operators using the AdS/CFT correspondence. These analyses allow us to picture 
the interacting, thermal field theory as a plasma medium. As discussed in §3.3, it is useful 
to ask how do probe particles introduced into such plasmas behave? This question is not 
only interesting from a theoretical viewpoint, but also from a pragmatic standpoint. For 
instance, in the case of the QGP, one is interested in knowing how much energy is lost by 
a quark produced in the deep interior of the plasma as it traverses outward. There, one 
does not have reliable computational methods owing to the strongly coupled nature of the 
plasma, but as we shall see, in the AdS/CFT framework one can again distill this question 
to a simple classical computation. In this section we will explore the various attempts in the 
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literature aimed at addressing this question, starting with the balhstic motion of quarks in 
the plasma and then turning to a discussion of the stochastic Brownian motion of stationary 
probes. 

4.3.1 Energy loss and radiation of moving projectiles 

To understand the energy loss of probe particles in plasma medium, we introduce an external 
probe in the form of an external quark or meson (for non-abelian plasmas) into the medium. 
Such probes are holographically modeled by an open string in the bulk geometry; here the 
geometry of interest is an asymptotically AdS black hole, which as we have described above 
provides the thermal medium. Let us understand the set-up in more detail. One of the open 
string end-points is pinned on the boundary of the AdS spacetime. Since this end-point 
carries the usual Chan-Paton index, it corresponds to the external quark we have introduced 
into the system. Heuristically, the external quark has a flux tube attached to it; in the 
holographic description one can view this flux tube as the open string world-sheet which 
extends into the bulk spacetime. Monitoring the motion of the external quark through the 
plasma thus amounts to studying the classical dynamics of string world-sheet in a black hole 
background. For mesonic probes we consider open strings with both end-points stuck on the 
boundary of AdS^+i. One could also consider other probes such as monopoles or baryons 
(which are heavy in the large- iV limit); these would correspond to D-branes living in the 
bulk. 

The first steps to understand the energy loss for probes in plasmas holographically were 
carried out in the seminal papers [87, 88, 89, 90, 91, 92, 93, 94], by considering the dynamics 
of probe strings as described above. A brief summary of these accounts can be found for 
instance in [9]. The general philosophy in these discussions was to use the probe dynamics 
to extract the rates of energy loss and transverse momentum broadening in the medium, 
which bear direct relevance to the physical problem of motion of quarks and mesons in the 
quark-gluon plasma. 

The simplest computation of the energy loss was originally considered in [87, 89]. The 
idea was to examine a probe in the boundary moving with a constant velocity v under the 
influence of an external force, applied so as to compensate for the frictional force acting on 
the probe and to maintain a steady state. From this steady state speed one can then recover 
the friction constant 7. On the bulk side, the problem reduces to a classical solution of 
Nambu-Goto action for the string, with one end-point being pulled with velocity v along the 
AdS boundary. The constancy of velocity at the boundary is again maintained by an external 
force, in this case a constant electric field that drags the string end-point. By examining 
the classical solution of the Nambu-Goto action with these boundary conditions, one finds 
that while the quark moves forward with constant velocity, the bulk string world-sheet trails 
behind, see Fig. 2. Since there is no natural place for the string world-sheet to end in the 
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Fig. 2: The trailing string solution in Schwarzschild-AdS4 spacetime. The curves are drawn for differing 
values of the quark velocity v (specifically from right to left, v = 0, 0.15, 0.5, 0.8, and 0.99999) 
on the boundary. 



bulk spacetime, it simply dips into the horizon. In particular, this implies that the classical 
world-sheet of the string itself has an induced horizon (which for non-zero velocity will always 
be outside the spacetime event horizon), a fact that will be important when we address the 
stochastic motion of the boundary end-point in §4.3.2. 

To see the construction in more detail, consider the planar Schwarzschild-AdSrf+i black 
hole (2.11); we are looking for a classical solution to the Nambu-Goto action 

'5ng = j c^Ve^ V-det(5(induced) , (4.7) 

with cr" being the world-sheet coordinates (t, a). The computation is easily carried out in 
static gauge t = t and a = r. In order to track the motion of a quark on the boundary, it is 
convenient to use the ansatz 

x\t,r) =vt + ^{r) , x' = 0\/i^l, (4.8) 

resulting in the Lagrangian density: 



£ = -^/l + rV(r)e(r)2--^ , (4.9) 

with /(r) given by (2.11). The solution is obtained straightforwardly by noting that the 
conjugate momentum to ^, tt^ = |^ is conserved and constant, which allows one to solve 

This is however not captured by the constant-t snapshot of the dragged string represented by Fig. 2 
where only the static region outside the horizon is visible, so the string looks like it 'freezes' onto the horizon, 
analogously to the "frozen star" picture of a collapsing black hole. 
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^(r) via quadratures 
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^'W = ^^W^77^^ 2- (4-10) 



In fact, requiring that the induced metric be timehke and non-degenerate fixes vr^ completely. 
Noting that < /(r) < 1 for r G [r+,oo) and f < 1 implies that the numerator of (4.10) 
can vanish somewhere in the bulk of the spacetime. At this point, it must be that the 
denominator also vanishes. This then constraints vr^ to take on the value: 

"^ = ^+(13^- (4.11) 

Given the classical string configuration obtained by solving (4.10), one can compute the 
rate of energy loss by looking at the momentum fiow along the string world-sheet. One 
simply has to compute the net fiux of this world-sheet momentum down the string. This 
results in [87, 89, 95] 



2 



AdS ^ _i 



mv 



P' = ^==^ ■ (4.12) 



dt d^a' (l-t;2)2/rf' ^ ^/i^ 

where we have translated the result in terms of the temperature T which is related to r+ 
via (2.14) and reinstated the AdS^+i length scale. In the non-relativistic limit, v <^ 1, this 
means that the friction constant is 

7o''="' = "^^l^^ . (4.13) 

a' m 

If we use the Sutherland-Einstein relation (3.21),^*^ we obtain the diffusion constant 

The calculation described above breaks down at very large velocities: physically, the force 
required to keep the quark moving sufficiently fast becomes so large that the quark anti-quark 
production is unsuppressed [94]. In [87] the more involved analysis of actually trying to see 
a moving quark slow down explicitly was also undertaken. This involves solving the time 
dependent equations arising from the Nambu-Goto action (4.7) with the initial conditions 
of non-zero velocity. 

In this context, there is an interesting puzzle: a-priori one would have expected that any 
particle moving through a plasma will lose energy to the medium and slow down. While we 
have seen that this occurs for quarks, described by the trailing string configuration (4.10), 



The physically correct sign of the square root corresponds to the bulk piece of the string world-sheet 
trailing the quark on the boundary. 
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Note that, as explained around (3.22), the relation (3.21) does not depend on d. 
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there is no mechanism for energy loss in the case of a meson moving through a plasma. A 
meson, modeled as a quark anti-quark pair corresponds to a configuration of an open string 
with both ends stuck to the boundary; in such situations one finds 'no-drag' solutions where 
the string simply dangles down into the bulk independently of the direction of v. More 
pertinently, for small separations between the quark and anti-quark (such that the meson 
remains in its bound state), the string stays above the horizon and as a result does not lose 
energy to the black hole. This feature was observed in different contexts in the analysis 
of [96, 97, 98]. Likewise it was also noticed in [99] that baryonic particle also propagate 
without being subject to drag (these authors also derived the drag force on k-quarks and 
gluons). These probes which seemingly do not suffer from dissipation are to our knowledge 
rather poorly understood. Analysis of the velocity dependence of the screening length in 
the plasma for mesons and baryons was explored in [97, 98, 93, 100]. Another mechanism 
for quark energy loss in the medium based on Cherenkov radiation was proposed recently in 

[101].22 

Note that deceleration of the quark provides yet another mechanism for its energy loss, 
induced by radiation due to the quark's deceleration. The interesting question of interplay 
between the two distinct energy loss mechanisms - namely medium-induced (i.e. drag) and 
acceleration-induced - has been explored in [102, 109], suggesting that these two effects may 
interfere destructively. More specifically, the authors consider a quark undergoing a circular 
motion at constant angular velocity. One advantage of such setup is that it can be treated as 
approximately stationary configuration while nevertheless incorporating acceleration (and in 
fact the classical world-sheet calculation for circular acceleration remains valid well into the 
acceleration-dominated regime), in contrast to the above-mentioned cases of linear motion. 
Interestingly, [109] find that depending on the angular velocity u and radius v/u of the 
quark's circular trajectory, the energy loss is dominated by either by drag force acting as 
though the quark were moving in a straight line, or by the radiation due to the circular 
motion as if in absence of any plasma, whichever effect is larger, with continuous crossover 
between these regimes occurring at w ~ TrT(l — f ^)^/^. 

Once one has an understanding of the basic mechanism for energy loss in the holo- 
graphic set-up, one can attempt to study the detailed response of the projectile in the plasma 
medium. One would naturally expect that a moving projectile produces a wake behind it; at 
large distances from the projectile one can employ hydrodynamics to study the behaviour of 

In fact, the gravitational dual is analogously puzzling. A dragged string above a static black hole is 
equivalent to a static string above a boosted black hole, since only the relative velocity matters. For the 
latter configuration, the boost of the black hole produces an ergoregion above the horizon, and it is easy to 
see that the string cannot dip into this ergoregion. However, one might naively expect the frame-dragging 
effect to extend past the ergosurface and affect the string, bending it in the direction of the boost. Yet, as 
manifest from the solutions, no such effect takes place! 

We should also remark that there have been attempts to study linearly accelerating particles and their 
associated radiation in the AdS/CFT context; see [102, 103, 104, 105, 106, 107, 108]. 
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this disturbance, but in general, high-frequency modes will be excited near the source of the 
perturbation. A convenient way to probe the physics of a moving projectile is to monitor 
the expectation value of the energy-momentum tensor: 

with projectile without projectile • 

(4.15) 

The projectile as before is modeled as an external quark in the field theory and corresponds 
to the motion of the string in the black hole background. Now in addition to the classical 
motion of the string we also want to compute the back-reaction on the geometry due to this 
external string and its influence on the boundary values of the metric (which in turn are 
related to the stress tensor expectation value). The first steps in this direction were taken in 
[110] where the expectation value of the operator dual to the dilaton was considered. This 
was then extended to the calculation of the stress tenor expectation value (4.15) in a series 
of works, [111, 112, 113, 114, 113, 115, 116, 117]. For a good review of this development we 
refer the reader to [118]. 

The essential idea here is to look at the perturbations about the trailing string configura- 
tion and ascertain from this the source of the bulk graviton. One then solves the (linearized) 
bulk Einstein's equations including the back-reaction of this source. Using the asymptotic 
behaviour of the solution one can finally read off the boundary stress tensor. This stress ten- 
sor exhibits many characteristic features that are expected for a projectile moving through a 
dissipative medium. For example, when the projectile moves faster than the speed of sound 
through the medium, it produces a Mach cone at a particular angle, leaving behind it a wake 
in the plasma (see [118] for further details). 

Using the setup of [109] describing quark in circular motion with constant angular ve- 
locity, [119] examine (part of) the corresponding boundary stress tensor induced by metric 
deformation due to the trailing string in the bulk. In particular, they compute the energy 
density and angular distribution of the power radiated by the quark. Unlike the previously- 
mentioned cases, they focus on the plasma medium at zero temperature, so that the bulk 
dual is described by a string trailing in pure AdS. A snapshot of a string is plotted in Fig. 3 
for various values of the quark velocity. The string solution looks like a rigidly-rotating 
spiral flaring out into the bulk, which in fact induces a horizon on the world-sheet despite no 
horizon being present in the underlying geometry (such horizon generation was previously 
discussed in [102] in more general context). ^'^ 

Rather curiously, [119] discover that this strong-coupling calculation gives very similar 
results to those at weak coupling, and closely resembles synchrotron radiation produced by an 
electron in circular motion in classical electrodynamics. Specifically, the radiation is emitted 
in narrow beam along its velocity vector, with opening angle ~ a/1 — w^. Surprisingly, 
despite strong coupling, the pulses of radiation propagate without broadening. While this 

The corresponding world-sheet temperature Tws is related to the Unruh temperature Tjj of the quark 
as Tws — Tu Vl — v^. 
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Fig. 3: The spiral string solution in AdSs spacetime corresponding to quark undergoing a circular 



motion. The curves are drawn for differing values of the quark velocity v on the boundary 
(values of v and color-coding same as in Fig. 2). The small black circle on the top is the 
quark's trajectory on the AdS boundary; the vertical direction corresponds to the bulk radial 
direction. 



is all explicitly derived from the bulk description, even on the gravity side it remains rather 
puzzling that the metric perturbation due to the trailing string remains so sharply localized. 
Indeed, this is rather counter to the naive UV/IR intuition, which would lead one to expect 
that the part of the string situated deep in the bulk should produce rather diffuse boundary 
stress tensor. 

The absence of broadening in the radiation pattern at zero temperature is in sharp 
contrast to the corresponding behaviour in the finite temperature case, where we expect 
that the beam of radiation propagating through thermal plasma slows down to the speed 
of sound and ultimately thermalizes. This type of behaviour was indeed observed in e.g. 
[120, 121] in the context of light quarks and mesons being released from rest. Here the 
authors start with initial state corresponding to a highly energetic quark anti-quark pair 
and follow its time evolution^'^ into two jets. Whereas at zero temperature the jets travel 
and spread forever, at finite temperature they stop within a finite distance Axmax, which 

scales with the cube root of the energy, Axmax ~ ( Tfryt ) 



A completementary approach to computing the energy lost by moving quarks in a non- 

Unlike heavy quarks, whose energy loss rate depends only on their velocity and the characteristics of 
the medium (i.e. the t Hooft coupling A and the temperature T of the plasma through which the quark is 
moving), light quarks are in addition influenced by details of their initial conditions, specified by the gauge 
field configuration. However, a more universal quantity, which is almost insensitive to the initial profile of 
the string, is the maximum distance Axmax traveled by a quark with initial energy E. 
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abelian plasma is to relate the energy loss to the so called "jet- quenching parameter" [122, 
123]. This parameter can be extracted from the expectation value of a light-like Wilson loop: 



( W{C) ) ~ exp 



4^2 



1 



qL- 



) 



(4.16) 



where the contour C is taken to be a rectangle of coordindate length L~ in the light-like 
direction and L in the transverse spatial directions with L" ^ L.^^ The computation of the 
expectation value of Wilson loops is well known from the early days of AdS/CFT [124, 125] 
and requires computing the area of a string world-sheet ending on an appropriate contour in 
the boundary of AdS. The first calculations of the jet-quenching parameter in A/" = 4 SYM 
were done in [88], where the authors found that 



There is by now a large literature exploring aspects of jet quenching in various set-ups and 
we refer the reader to the original papers for more details. 

The discussions so far have focussed on the dissipative aspects of the non-abelian plasma. 
As we have seen the AdS / CFT correspondence allows one to extract the basic properties of 
such media; specifically, by examining the classical energy loss by probes into black holes 
one recovers the frictional characteristics of such plasmas. 

4.3.2 Brownian motion in AdS/CFT 

Thus far we have focused on the classical motion of a probe particle in a thermal plasma 
medium. As expected, the particle loses its kinetic energy to the plasma and slows down. In 
the dual gravitational description of this frictional motion, this is mimicked by energy loss 
into the black hole horizon. However, given that we have a particle in a thermal medium, 
one should expect that it not only slows down, but also undergoes random motion due to the 
thermal fluctuations of the plasma. This is understood as Brownian motion and, as we have 
explained in §3.3, can be modeled in terms of a Langevin equation, which arises naturally in 
the context of kinetic theory. It is therefore inviting to ask whether these random fluctuations 
can be captured in the gravitation description, and if so, what is their physical origin? 

This question is not only interesting for academic reasons, but as described above, the 
holographic models provide useful toy examples to study the physics of the QGP formed 
in heavy-ion collisions. In that context, the motion of an external quark in the QGP is 
assumed to be described by a relativistic Langevin equation [126]. As reviewed in §3.3, in its 
most basic form the Langevin equation is parametrized by two constants: the friction (drag 
force) coefficient 7 and the magnitude of the random force k. However, in the relativistic 

^'^Note that the physical proper length along a light-like direction is of course always zero. 



9A^=4SYM — 




(4.17) 
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case, the random force has different magnitudes kt and kl in the directions transverse and 
longitudinal to the quark momentum p, which in the non-relativistic limit p — > become 
equal, kl = kt- The parameters 7 and kl are related to each other by the Einstein relation, 
under the assumption that the Langevin dynamics holds and gives the Jiittner distribution 
g-/3E ^YiQ other hand, kt is an independent parameter [126, 92]. Furthermore, random 
force is usually assumed to be white noise, i.e, with a delta function auto-correlation. 

One would therefore like to understand to what extent this picture works in the holo- 
graphic setting. In early works on the subject, [91, 92, 94] computed the random force 
strength k^- The calculations of [91, 94] were carried out by expressing the change in the 
density matrix associated with the probe (and therefore k) in terms of a Wilson loop average. 
This Wilson loop was taken to lie along a Schwinger-Keldysh contour and corresponds to a 
source of the random force operator. On the other hand, [92] computes the random force 
strength by looking for fluctuations around the trailing string solution discussed in §4.3.1. 

In [95, 127] the complete story for the holographic Brownian motion was worked out in 
detail. This was further generalized in [128] to the relativistic Langevin equation. As we have 
alluded to earlier, classical string solutions in the black hole backgrounds have an induced 
world-sheet metric which has a horizon. In general, a world-sheet horizon could correspond 
to any of the following: 

• to the bulk spacetime horizon, which is the case for static strings in static black hole 
spacetimes 

• to a bulk ergosurface, which occurs for a stationary string solution like the trailing 
string or in the case of a stationary black hole; the ergosurface occurs at the location 
where the asymptotic timelike Killing field has a vanishing norm, i.e. Qu = 0, 

• or simply when the string is forced to move too fast, which occurs for instance when 
it flares out while rotating at constant angular velocity as in [119]. 

The key physical result of [95] was that the fluctuations of the classical string due to the 
curved spacetime particle production, i.e.. Hawking radiation, associated with the world- 
sheet horizon was responsible for the random force on the probe particle. In fact, this 
connection was also noted in some of the earlier works on the subject, notably [92] and 
[129]. Explicit expressions of this were worked on by assuming a thermal spectrum of the 
Hawking quanta in the case of BTZ black hole in [95], which had the advantage that the mode 
functions could be analytically obtained. Higher dimensional examples were also discussed in 
[95] and [127], the latter working directly with the maximal analytic extension of an eternal 
Schwarzschild-AdSrf+i black hole. The extension of these results to the trailing string solution 

Implicit in these statements is an assumption that one is working in static gauge, where the world-sheet 
time coordinate is identified with the timelike Killing field of the stationary asymptotically AdS^+i black 
hole. 
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was described in [128, 130]. In addition [131] examined the motion of a heavy quark in a 
toy model of a dynamical black hole, the so called conformal soliton geometry,^'' which we 
will revisit in §6.1. Similar studies for the case of an accelerating quark in the vacuum 
was recently considered in [107]. Futhermore, in the authors of [133, 134] have examined the 
details of Langevin dynamics under the improved holographic QCD framework; in particular, 
these works compute the full Langevin correlator in non-conformal models, with the aim of 
using these results for studies of QGP. An important caveat to bear in mind for applications 
to QGP is that one generically has to consider the entire correlator and can not simplify the 
dynamics to the basic Langevin dynamics discussed in §3.3. 

The simplest setting to study the stochastic motion is to monitor the fluctuations about 
the static straight string. For simplicity let us consider the planar AdS^+i black hole (2.11) 
and take the string world-volume to be = 0. Expanding the fluctuations of the string to 
quadratic order about this solution, we obtain from the Nambu-Goto action (4.7): 

1 /" ,9 , — ,„/ dx^ dx^ 
where the induced metric ■jf^^ in static gauge is 

dsLuce, = -r' f{r) de + , /(r) = 1 - -± (4.19) 

and the kinetic terms of the scalars are determined by gij = 6ij. The equation of motion 
for the transverse scalars on the probe string world-sheet then takes the form: 

- d^x' + f{r) dr [r^ f{r) drx') = (4.20) 

This wave-equation for the non-minimally coupled transverse scalar x{t,r) is solved by 
standard mode expansion x^{t, r) = e^**^* ^*(?"); in general d this is somewhat complicated by 
the fact that the explicit mode solutions are not known. However, for d = 3, i.e., for BTZ 
black holes one can solve this wave equation explicitly. Analysis of the wave equation in the 
near-horizon region r ^ r+ reveals that the linearly independent solutions can be taken to 
be the ingoing and outgoing waves f(±)(r). Thus one can generally write the solutions after 
implementing the appropriate boundary conditions on the horizon and at infinity in terms 
of a mode expansion: 

x\t,r) = J2 h>L(i,r) + (aL)^«L(i,r)*] 

ui>0 

ul{t,r) = A{e^+\r) + BC^'\r)) e-''^' (4.21) 

The conformal soliton spacetime was first described in [f32] and is simply the global Schwarzschild- 
AdSd+i spacetime considered in a single Poincare patch. From the field theory point of view, it corresponds 
to a conformal mapping of a thermal stress tensor on S'^""'^ x R to R''"^'^. We will have occasion to discuss 
this spacetime in the context of the fluid/gravity correspondence in §6.f . 
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with the creation and annihilation operators and aj^ satisfying the standard commutation 
relations: 

k, <'] = = 0, k, {ai,y] = 5^^, 5^^ (4.22) 

So far we have outlined how to write down the canonical quantized description of the 
fluctuations of the transverse excitations of the static string solution. The essential feature 
we now wish to exploit is that the asymptotic behaviour of the fluctuations (considering a 
regulated AdS spacetime with cut-off at r = 

xUt)= lim x'{t,r) , (4.23) 

captures the motion of the probe quark in the boundary plasma. Monitoring the dynamics of 
^oo(^) allows one to compute the details of the stochastic motion of the probe. In particular, 
one has a precise mapping between the correlation function of the quark in the boundary 
plasma and the correlators of the quantum operators a^, i.e., 

(ti)xoo(t2)...) (4.24) 

This mapping in particular implies that, were we able to make a very precise measurement 
of the correlators of the end-point of the string, then one could in principle make a precise 
measurement of the state of the radiation emanating from the black hole (for after all the 
correlation function in the bulk depends on the state of the theory we are considering). 

However, if we are interested in understanding the stochastic process at a semi-classical 
level, all one needs to do is to evaluate the bulk correlation function in the natural Hartle- 
Hawking state for the fluctuating flelds x^{r,t). It is a well known fact that the Hartle- 
Hawking state corresponds to thermal equilibrium and it therefore follows that the outgoing 
mode correlators are determined by the thermal density matrix 

p = e-^\ H = J2^Ky< (4.25) 

aj>0 

Using this one can compute the mean squared displacement (: [xoo(t) — Xoo(O)]^ :) (for 
simplicity monitoring only one of the scalar flelds x\ say x^ = x). This allows one to extract 
the random force correlation function k,{u). Similarly, one can compute the admittance 
p{u}) deflned in (3.29) by examining the behaviour of the string end-point when subject to 
an external force. As one expects the results conflrm the fluctuation-dissipation theorem. 
We refer the reader to [95] for the precise expressions. It is more useful to record here the 
various time-scales obtained from the holographic computation: 

a' m 1 1 

^relax ~ ^2 ^ ' ^^oll ~ 7f^ , ^mfp 7=— (4.26) 

which in particular implies that 

^relax a' 171 171 tcoU 



?2 



^coll -^AdS^ vAT ^mfp 



A > 1 (4.27) 
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which clearly demonstrates the deviation from the hierarchy mentioned in §3.3. One can 
understand this as a clear signal of the strongly coupled nature of the plasma dual to the 
black hole geometry, which obviates the wide separation of time-scales pertaining to kinetic 
theory. 

In [127] the stochastic motion described above was derived using an extension of the 
Schwinger-Keldysh formalism. Instead of working with the static region of the Schwarzschild- 
AdS black hole, one considers maximally extended spacetime in Kruskal coordinates, and 
utilizes the Lorentzian AdS/CFT prescription of [57] to obtain the retarded Green's func- 
tions. For instance the random force correlator K,{t) is given in terms of the retarded Green's 
function: 



This calculation leads to pretty much the same physical picture of the random motion as the 
one described above. There are some essential differences in the derivation, since the real 
time formalism employed in [127] is intrinsically tied to the Lorentzian geometry, whereas 
the Hartle-Hawking state used in the derivation of [95] is essentially an equilibrium thermal 
computation (which therefore can be carried out in the Euclidean black hole geometry). 

The description of the stochastic process described so far corresponds to the behaviour of 
the quark end-point on the boundary of the AdS spacetime (rather on a regulated boundary). 
An interesting question is to ask whether one can use this picture to further the membrane 
paradigm picture of the black hole. The main philosophy of the membrane paradigm [135] is 
that, as far as an observer staying outside a black hole horizon is concerned, physics can be 
effectively described by assuming that the objects outside the horizon are interacting with 
an imaginary membrane, which is endowed with physical properties, such as temperature 
and resistance, and is sitting just outside the mathematical horizon. 

In [95, 127] the stochastic properties of this fiducial membrane were explored, mainly by 
'integrating out' the string world-sheet all the way down to a stretched horizon located at 
r = + e. It was found that the heavy quark diffusion constant on this membrane was 
exactly the same as the asymptotic value (4.14). This is in fact similar to the results on 
the hydrodynamic properties of the membrane (which we will revisit in §5) as discussed for 
instance in [59]. 

A crucial question about the membrane paradigm is to understand the microscopic struc- 
ture of this stretched horizon in the context of string theory. In [136, 137], Susskind and 
collaborators put forward a provocative conjecture that a black hole is made of a fundamen- 
tal string covering the entire horizon. Although this picture must be somewhat modified 
[138] since we now know that branes are essential ingredients of string theory, it is still an 
attractive idea that, in the near horizon region where the local temperature becomes string 
scale, a stringy "soup" or "cloud" of strings and branes is floating around, covering the entire 
horizon. Ref [95] and more recently [139] have attempted to find precise characteristics of 





(4.28) 
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this membrane by trying to use the characteristics of the stochastic motion of strings ending 
on the stretched horizon. In particular, one could imagine the string IR end-point interact- 
ing with the membrane halo of the black hole. Any such model requires that we be able 
to match the mean free time between the interactions of the string end-point with the halo 
with the results derived above for the stochastic motion. The calculations in [95, 139] seem 
to suggest a non-trivial dependence of the interaction probability on Qs and Is- 

5 Fluid/gravity correspondence 

Thus far our understanding of physics out of equilibrium has been restricted to the linear 
response regime. In the holographic set-up this regime is captured by the study of linearized 
fluctuations about AdS black hole backgrounds. We have in particular seen that the retarded 
correlators of local gauge invariant operators can be efficiently computed in the classical 
gravity approximation and the results agree with general expectations from the field theory. 

A particular application of the linear response story is the use of Kubo formulae to extract 
transport coefficients of strongly coupled gauge theory plasmas. As described earlier, these 
studies led to the fascinating bound on the ratio rj/ s. A natural question is whether one can 
derive non-linear hydrodynamics from the holographic set-up, by relaxing the constraints 
of linearized fluctuations. A-priori this sounds hard, for the problem surely translates to 
understanding non-linear fluctuations around a black hole background. Nevertheless, as 
we now shall describe, this is indeed tractable and leads us to a natural relation between 
Einstein's equations and a relativistic generalization of the Navier-Stokes equations, which 
has come to be known as the fluid/ gravity correspondence, [22]. 

The key feature we will exploit in making this connection is the general idea that the 
systems in local thermal equilibrium should in a suitable infra-red limit admit a hydrody- 
namic description. As familiar from everyday experience, this framework can accommodate 
systems with large time-dependence, as long as they are in local thermodynamic equilibrium; 
proximity to global thermodynamic equilibrium is no longer required. 

More speciflcally, fluid dynamics is the continuum effective description of any (interact- 
ing) microscopic quantum fleld theory. In order to meaningfully describe the system in terms 
of the fluid variables, the fluid description requires that the system achieves local thermo- 
dynamic equilibrium. This means that the regime of validity where such a description is 
valid requires that the scale of variation of the dynamical degrees of freedom, L, be much 
larger than the microscopic scale £mfp, typically set by the temperature, T (or the local en- 
ergy density). In this long-wavelength approximation, local equilibrium then demands that 
LT = i > 1. 

To keep the discussion and formulas as clean as possible, we will restrict attention to 
the simplest case of uncharged conformal fluid on four-dimensional Minkowski spacetime 
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R^'^ (though when sufficiently compact, we will quote the d dimensional results). Several 
generalizations to this setup will be mentioned in §5.4. For basic details on fluid dynamics we 
refer the reader to [141]. Relativistic fluids are described in [142], aspects of the fluid/gravity 
correspondence are reviewed in [24, 25] and applications of relativistic fluids to heavy-ion 
collisions in [143]. 

Before proceeding, we make few remarks on notation: the bulk metric will be denoted 
by Qmn with the capital Latin indices taking values over the d + 1 bulk dimensions; we will 
separate the coordinates into the radial coordinate r and the remaining 'boundary coordi- 
nates' x^, where the /i index ranges over the d boundary directions (which includes time). 
The stress tensor in the boundary theory is denoted T^'^ , and in writing its conservation 
{y ^T^^^ = 0), the is the covariant derivative with respect to the boundary metric r/^;/.^^ 

5.1 Background 

We begin with a brief review of conformal fluid dynamics, proceed to discuss the dual 
gravitational solutions, and then motivate the construction of the explicit mapping between 
them. 



5.1.1 Conformal fluid dynamics 

A conformal fluid is characterized by a traceless symmetric stress tensor, which in d spacetime 
dimensions has ^^i^^til — i degrees of freedom, along with a collection of charge currents (which 
for simplicity we have set to zero). In a fluid dynamical characterization of the same system, 
the number of basic degrees of freedom is drastically reduced. The conformal invariance 
flxes the equation of state, thereby determining the pressure in terms of the energy density, 
which can in turn be expressed in terms of the temperature. Hence the basic variables are 
the local temperature T{x) and velocity u^{x) (unit-normalized so that u^u^ = —1) , which 
constitute just d degrees of freedom. 

The equations of fluid dynamics are then simply the equations of local conservation of 
the stress tensor (as well as the charge currents in more general situations), supplemented 
by constitutive relations that express these currents as functions of the fluid dynamical 
variables. As fluid dynamics is a long wavelength effective theory, such constitutive relations 
are usually specifled in a derivative expansion, like in any effective fleld theory. At any given 

^® The most comprehensive discussion of conformal fluids in arbitrary dimensions can be found in [140]. 
Oftentimes the reader will encounter an ordinary partial derivative 9^ in the following: this is because 
in Cartesian coordinates of flat spacetime, V^i = d^. However, this is convenient more generally as well: 
since the long-wavelength approximation engenders an ultra-locality in the equations one solves (along the 
boundary directions), wc can use boundary derivatives as just ordinary partials; however, our expressions 
can of course be easily 'covariantized' by replacing 9^ by V^. We hope this does not cause much confusion. 
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order, thermodynamics plus symmetries determine the form of this expansion up to a finite 
number of undetermined coefficients. In general, the coefficients can be obtained either from 
measurements or from microscopic computations. However, as we will see, in the present 
framework these coefficients are fully determined by the gravity side (which in a sense knows 
about the microscopies of the boundary field theory). 

Purely based on the symmetries, we can then write down an expression for the stress 
tensor of a (i-dimensional conformal fluid, which is a local functional of the temperature and 
velocity flelds: 



The flrst two terms describe the ideal conformal fluid stress tensor, while vr^^^jp^^j^^ incorpo- 
rates all the dissipative terms, a here sets the overall normalization of the stress tensor. As 
variations of T{x) and u^{x) are small, we can expand tx^" in a derivative expansion c?^ = ^ 
in the boundary directions; the leading term will turn out to be proportional to the shear 
viscosity. The dynamical content of the fluid equations is encoded in the conservation of the 
stress tensor 



Fluid dynamics viewed in this derivative expansion constructs an effective fleld theory for 
the slowly varying modes T[x) and analogously to the chiral Lagrangian for pions. 

In general one should write down all possible terms at a given order in the derivative 
expansion consistent with the symmetries of the problem (and modulo lower order conserva- 
tion equations). There are scalar functions of the thermodynamic variables multiplying these 
operators: these are the transport coefficients. For instance, we have indicated in (5.1) the 
explicit terms occurring at zeroth order (the ideal fluid terms). At the next order, there is one 
term for a conformal fluid (proportional to the shear viscosity), while a non-conformal fluid 
would have two. At higher orders we get more terms; a discussion of terms allowed at second 
order was first undertaken in [83] in four dimensions and higher dimensional conformal fluids 
were discussed in [144, 140] and a nice account of non-conformal fluids can be found in [145]. 
We should note that conformal fluids are especially simple, as Weyl covariance provides a 
very useful tool to classify the possible terms in the fluid stress tensor. Using a manifestly 
Weyl-covariant formalism [146], many of the expressions derived in [22] simplify. We refer 
the reader to [140] for compact expressions in diverse dimensions. 

5.1.2 Gravity in the bulk 

We now turn to the gravitational solutions in asymptotically AdS spacetime. Motivated by 
the AdS/CFT correspondence, we will consider two-derivative theories of gravity with an 
AdSrf+i "vacuum", such as the IIB SUGRA on AdSs x S^. As mentioned in §2.1, the solution 
space has a universal sub-sector, pure gravity with negative cosmo logical constant, for which 





V^T^"^ = . 
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the bulk field equations are simply Einstein's equations, 

Emn = Rmn — - Rqmn + ^dMN = . (5.3) 

(Note that taking -RAds = 1 sets A = —6 in five dimensions.) We will focus on this sub-sector 
in the long-wavelength limit. Apart from the pure AdSs solution, there is a 4-parameter 
family of solutions representing asymptotically- AdSs boosted planar black holes. We will use 
these solutions to construct general dynamical spacetimes characterized by fiuid-dynamical 
configurations. 

Roughly speaking, the fluid/gravity construction may be thought of as a "collective 
coordinate method" for black hole horizons. Recall that the isometry group of AdSs is 
5*0(4,2). The Poincare algebra plus dilatations form a distinguished sub-algebra of this 
group (one that preserves the boundary). Out of these, the 5*0(3) rotations and translations 
in world-volume R^'^ leave the static planar AdS black hole invariant, but the remaining 
symmetry generators, dilatations and boosts, act nontrivially on this solution, generating 
a 4-parameter family of boosted planar black holes, parameterized by the temperature T 
and the velocity u'^ of the brane. The construction effectively promotes these parameters 
to Goldstone fields (or collective coordinate fields) T{x) and u^{x), and determines their 
dynamics, order by order in the boundary derivative expansion. Note that this is distinct 
from linearization: we make no assumptions about the amplitudes of these slow variations. 



5.1.3 The fluid/gravity map 

Before proceeding to sketch the construction in more detail, we pause to stress an important 
point in mapping these long-wavelength gravity solutions to corresponding fiuid configura- 
tions. A well-known procedure of holographic renormalization (see e.g. [44] and [26]) links 
the boundary stress tensor to the behaviour of the bulk metric near the AdS boundary. 
Given any asymptotically AdS spacetime, we can read-off the induced stress tensor on the 
boundary, since the latter is related to the normalizable modes of the gravitational field 
in AdS. In particular, expanding the bulk metric in the Fefferman- Graham form near the 
boundary 2; = 0, 

,2 dz'^ + iri^^ + az'^ T„^) dw^" dw" 

ds = ^ — ^ , 

z^ 

the stress tensor is simply given by T^,y. Conversely, given a boundary stress tensor, there is 
a procedure to holographically reconstruct the bulk metric in a radial expansion around the 
boundary. 

Naively, this might seem puzzling: as mentioned above, a conformally invariant stress 
tensor in d dimensions has ^^^tii — 1 degrees of freedom. If any such stress tensor yielded a 
regular bulk spacetime, we would have a discrepancy between the fiuid side which has only 
d degrees of freedom and whose dynamics is correspondingly specified by only d equations. 
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and the gravity side that would seemingly allow more degrees of freedom. In other words, 
passing from a generic quantum conformal field theory stress tensor to the stress tensor of its 
effective description in terms of fluid dynamics constitutes a drastic reduction in the number 
of degrees of freedom required to specify the spacetime. How is this manifested in the bulk? 
The answer lies in regularity. As a series expansion around the boundary, the holographic 
reconstruction cannot guarantee that the metric does not become nakedly singular at some 
finite radial value in the bulk. In fact, for a generic stress tensor it will. The fluid/gravity 
construction demonstrates that the regular solutions are given precisely by such stress tensors 
which are fluid dynamical. Moreover, we claim that the gravity solutions thus constructed 
are the most general regular long-wavelength^° solutions to Einsteins equations with negative 
cosmological constant. They typically correspond to deformed and dynamical black holes; 
i.e. the solutions admit a regular event horizon which shields a curvature singularity. 

The heuristic picture of a generic evolution, on the two sides of the fluid/gravity corre- 
spondence, is as follows. Suppose we start with some generic high energy initial conditions. 
On the CFT side, the system quickly settles down to local thermodynamic equilibrium, 
whose bulk dual is described by a dynamical, non-uniform (planar) black hole. On both 
sides, such configuration is described by local velocity and temperature fields which exhibit 
slow variation in the boundary directions. The subsequent evolution is described by equa- 
tions of fluid dynamics on the boundary, which originate from Einstein's equations governing 
the bulk evolution. Finally, at late times, the system relaxes to global thermal equilibrium, 
given by a stationary state parameterized by a constant temperature and velocity. In the 
bulk, this is one of the well-known stationary solutions describing a planar black hole in AdS 
mentioned in the previous subsection and given explicitly below. 

The fluid/gravity construction specifically utilizes the fact that fluid dynamics is a long- 
wavelength effective theory. One writes Einstein's equations as a perturbative expansion 
in boundary derivatives (however keeping the exact radial dependence) to emphasize the 
expansion at small momenta. This allows us to solve the equations order by order in this 
boundary derivative expansion. In turns out that Einstein's equations at a given order 
implement the fluid stress tensor conservation equations at lower order. Therefore, order by 
order, we can use the lower-order fluid dynamical solution to construct the bulk metric, and 
then read off the corrected fluid dynamical stress tensor. In [22] , the boundary stress tensor 
T^^" and corresponding bulk metric qmn were constructed to second order in the boundary 
derivative expansion. This yields a map between fluid dynamics and gravity, which we now 
proceed to sketch in more detail. 

'^'^ Note, however, that there are regular solutions which do not fall into the long-wavelength category, such 
as small black holes in AdS, which correspondingly are not described by fluid configurations. 



47 



5.2 Iterative construction of bulk metric and boundary stress tensor 

The iterative procedure starts with the zeroth order configuration, corresponding to the 
global equilibrium, given in (2.13). As will shortly become clear we will implement a per- 
turbation expansion about this solution. In order for perturbation theory to make sense, it 
is important that we start with a seed metric in manifestly regular coordinates. 

A general fluid conflguration in local, but not global, equilibrium can be described by 
promoting the 4 parameters T and to physical fields dependent on the boundary coor- 
dinates x^, i.e., to T{x) and u^{x). If these fields vary slowly compared to the microscopic 
scale ^mfp, i-e. if 

for small e, the fiuid configuration still satisfies the conditions of local equilibrium. In each 
local domain of slow variation, which we refer to as tube, the bulk gravitational solution 
is approximately that of a uniform black brane. Remarkably, the bulk solution can be 
constructed by patching together these tubular domains! Of course, if we just replace 
and T in the metric (2.13) by T{x) and u^{x), the resulting metric (call it q^^n) ^^^^ 
longer solve Einstein's equations (5.3). Instead, the metric gfj^^ will need to be corrected 
by higher-order piece {q^mni etc.), which we can obtain iteratively as an expansion in e. We 
will find that the resulting corrected metric can be constructed systematically to any desired 
order, and is valid well inside the event horizon, thus allowing verification of its regularity. 
It is worthwhile to stress that the success of such a procedure rests on the fact the our seed 
metric g^^^ is manifestly regular on the horizon, since otherwise the expansion would break 
down near the coordinate singularity at horizon. 

To implement the construction algebraically, we express the line element in a boundary 
derivative expansion of the fields u^{x) and T[x), and use e as a book-keeping parameter 
(counting the number of x^ derivatives): 

oo oo oo 

9MN = E , T = J2e' T(^) , u, = J2e' . (5.4) 

fc=0 fc=0 fc=0 

The term gf]^ corrects the metric at the k^^ order, such that Einstein's equations will be 
satisfied to 0{e^) provided the functions T{x) and u^{x) obey a certain set of equations of 
motion, which turn out to be precisely the stress tensor conservation equations of boundary 
fiuid dynamics at C(e^~^). 

Specifically, we can obtain the equations for g^^^^ by substituting the expansion (5.4) into 
Einstein's equations (5.3), and extracting the coefficient of 0{t^). Schematically, these take 
the form 

'H[g^'\uf,T^%g'^'\xn = su (5.5) 

where 7^ is a second-order linear differential operator in the variable r alone and Sk are 
regular source terms which are built out of (yf*-"^ with n < k — 1. Since g'^^\x^) is already of 
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0{e''), and since every boundary derivative appears with an additional power of e, "H is an 
ultralocal operator in the field theory directions. Moreover, at a given x^, the precise form 
of this operator H depends only on the local values of T and but not on their derivatives 
at x^. Furthermore, the operator H is independent of k; we have the same homogeneous 
operator at every order in perturbation theory. This allows us to find an explicit solution of 
(5.5) systematically at any order. The source term Sk however gets more complicated with 
each order, and refiects the nonlinear nature of the theory. 

Bit more explicitly, the equations of motion split up into two kinds: Constraint equations, 
Er^i = which implement stress-tensor conservation (at one lower order), and Dynamical 
equations E^,^ = and Err = which allow determination of g^''\ We solve the dynamical 
equations 

g^''^ = particular(sfc) + homogeneous (H) 

subject to regularity in the interior and asymptotically AdS boundary conditions. Using the 
rotational symmetry group of the seed solution (2.13) it turns out to be possible to make a 
judicious choice of variables such that the operator "H is converted into a decoupled system of 
first order differential operators. It is then simple to solve the equation (5.5) for an arbitrary 
source Sk by direct integration. For the details of the procedure, as well as discussion of 
convenient gauge choice, etc., we refer the reader to the original work [22] or the review [24]. 

Instead, here we simply quote the result for the bulk metric and boundary stress tensor, 
corrected to first order in e. To first order the bulk metric takes the form 



ds'^ = —2u^dx'^dr + r'^{ri^y + [\ — f{r/'nT)]u^Uy)dx^dx'' 



+ 2r 



r 11 
— F{r/ TxT) + - u^Uy dxu^ - - u^dx (u^u^) 



dx'^dx", (5.6) 



where T{x) and u^{x) are any slowly- varying functions which satisfy the conservation equa- 
tion (5.2) for the zeroth order ideal fiuid stress tensor (2.17), the function F{r) is given 
by 



Fir) = r dx ^' + ^ + 1 1 



In I — 2 arctan(r) -|- tt 



x(x+l)(x2 + l) 4 V 
and a'^" is the transverse traceless symmetric part of d'^u'^ called shear, i.e. 

Note that the first line of (5.6) is simply the zeroth order (boundary-derivative-free) solution 
(2.13), whereas each of the terms in the second line have exactly one boundary derivative. '^^ 



Note that (5.6) does not have any terms appearing exphcitly, since by implementing the zeroth 
order stress tensor conservation, we have expressed the temperature derivatives in terms of the velocity 
derivatives. 
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The induced fluid stress tensor on the boundary, which can be easily obtained^^ from the 
bulk metric (5.6), is given by 

rj^i^v ^ ^4 y4 ^4 _^ ^i.u^ -2n^T^ a^"". (5.7) 

Here the first two (derivative-free) terms describe a perfect fluid with pressure (or negative 
free energy density) vr^T^, and correspondingly (using thermodjTiamics) entropy density 
s = An'^T^. The shear viscosity r] of this fluid may be read off from the coefficient of a'^" 
and is given by vr^T^. Notice that rj/s = 1/{Att), in agreement with the well-known result 
of [147]. 

5.3 Solution at second order 

In the previous section we have illustrated at first order how the iterative procedure can be 
implemented (in principle systematically to any order) to construct a generic long-wavelength 
solution. Such a procedure was carried out in [22], where the bulk metric and boundary stress 
tensor was calculated explicitly to second order in the boundary derivative expansion. In 
this section we will discuss the new physics which can be extracted from such a construction. 

Note that already at first order, the bulk metric (5.6) was a much lengthier expression 
than the boundary stress tensor (5.7). This remains true in general; in fact, already at 
second order the expression for the metric is far too unwieldy to write down here. In the 
following subsections, we will therefore only write the second order boundary stress tensor 
explicitly but indicate the bulk metric only schematically. 

5.3.1 The 4-dimensional conformal fluid from AdSs 

The second order stress tensor obtained from the gravity analysis is best expressed in terms 
of Weyl invariant operators. To do so it is useful to classify all the operators which are Weyl 
invariant at various orders in the derivative expansion. At first two orders in derivatives the 
set of symmetric traceless tensors which transform homogeneously under Weyl rescalings are 

■^^ One can compute the stress tensor using the general formula (2.16), which can be reduced in the current 
case to the simpler formula: 

Ti: = -2 lim (X^ - <5t') . 
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given to be^^ 

First order : a^'^ 
Second order : 1^''' = 2 Vacr^^, , Tg'' = C^aup u'^ 



where we have introduced a notation for the second derivative operators which will be useful 
to write compact expressions for the stress tensor below. The quantities involved in the op- 
erators above are constructed from the velocity derivatives, such as the acceleration a^, shear 
a^^ , etc.. These quantities are defined using the decomposition of the 4-velocity gradient 
V^u^ into transverse, traceless and trace parts, 

1 
3 

where expansion, acceleration, and vorticity, are respectively defined as: 

e = V^M^ , = u'^V^u^ , uj^f" = P^" P^^ V[«M/3] • 

In addition in four spacetime dimensions we also have the 'curl' of the velocity field 

p = e"^^^ V/3M^, (5.8) 

Armed with this data we can immediately write down the general contribution to the 
stress tensor as: 

nf,") = -2r/(T^" 

nf^") = r^r/Tf + z^Tr + AiTf + A2Tf + A3Tr . (5.9) 

There are therefore have six transport coefficients rj, r^r, n, Aj for i = {1,2,3}, which char- 
acterize the flow of a non-linear viscous fluid. 

For a fluid with holographic dual using the result of the gravitation solution (the pro- 
cedure to construct which is described in §5.2) one finds explicit values for the transport 
coefficients. In particular, for the A/" = 4 SYM fluid one has^^ [83, 22] 

(vrT) 



8tt s An 

2-ln2 _ 7] 

2tiT ' ^~ ^ 



5 



(5.10) 



\ - r7ln2 

■^■^ Various papers in the literature seem to use slightly different conventions for the normalization of the 
operators. We will for convenience present the results in the normalizations used initially in [83]. This is 
the source of the factors of 2 appearing in the definition of the tensors Ti, see [22] for a discussion. The 
derivative is a Weyl covariant derivative introduced in [146] which makes it easy to keep track of conformal 
transformation properties of various tensors. 

■^"^The result for generic M — 1 superconformal field theories which are dual to gravity on AdS5 xA's are 
given by simply replacing by the corresponding central charge of the SCFT. 
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where in the first hne we have used the standard entropy density for thermal Af = 4 SYM 
at strong coupling s = ^ N"^ to exhibit the famous ratio of shear viscosity to entropy 
density [15]. 

So far we have only discussed the behaviour of four dimensional relativistic fluids which 
have bulk duals as asymptotically AdSs black hole spacetimes. The story is readily general- 
ized to other dimensions d > 3. The general analysis in AdS^+i of carried out in [144, 140] 
in fact allows us to write down the transport coefficients described above in arbitrary dimen- 
sions in nice closed form.^^ The transport coefficients for conformal fluids in li-dimensional 
boundary Bd are: 



1 fAirV'^ r] 1 



d 



471 T 



1.. ^2 
1 + - Harmonic - — 1 
d \ d 



d rj 
^ ~ 27r {d - 2) T 



(5.11) 



d 7] ,1 (2 \ 7] 

Ai = --, = -Harmonic (--1-, A3 = 



where Harmonic fx) is the harmonic number function 
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5.3.2 The spacetime geometry dual to fluids 



Let us now turn to discuss the bulk geometry obtained at the second order in boundary 
derivative expansion (whose first order part is given by (5.6)). As mentioned previously, this 
bulk solution is 'tubewise' approximated by a planar black hole. This means that in each 
tube, defined by a small neighborhood of given x^, but fully extended in the radial direction 
r, the radial dependence of the metric is approximately that of a boosted planar black hole at 
some temperature T and horizon velocity u^. These parameters vary from one position to 
another in a manner consistent with fiuid dynamics. Our choice of coordinates is such that 
each tube extends along an ingoing radial null geodesic; see Fig. 4. Apart from technical 
advantages, this is conceptually rather pleasing, since it suggests a mapping between the 
boundary and the bulk which is natural from causality considerations. 

It is worth stressing that although we refer to the metric written in (5.6) and its second 
order extension presented in [22] as "a solution" in the singular form rather than the plural, 
these expressions actually correspond to not just a single solution or even a finite family of 

■^^For d — i the general results were initially derived in [148]. See also [85] for determination of some of 
these coefficients in d = 3 and d = 6. 

■^^The harmonic number function may be in fact be re-expressed in terms of the digamma function, or 
more simply as 

■ r ^ r'fx) 

Harmonic(a:) = 7e + -^^^^ , 

where 7e is Euler's constant. 
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Fig. 4: The causal structure of the spacetimes dual to fluid mechanics illustrating the tube structure. 

The dashed line denotes the future event horizon T-L'^ generated by while the shaded tube 
indicates the region of spacetime over which the solution is well approximated by a tube of 
the uniform black brane. 

solutions, but ratlier a continuously infinite family of solutions, specified by the four functions 
T(x) and u^{x) of four variables. The flip side of the coin is that while very general, such 
a metric is not fully explicit: in order to be so, we need to use a given solution to fluid 
dynamics as input. 

However, even in the absence of the explicit functional dependence of T{x) and u'^{x), 
it is possible to extract certain salient features of any such geometry. The most important 
feature of our geometry is the presence of an event horizon. In [149] we have demonstrated 
explicitly that the event horizon is regular, and determined its location in terms of the 
functions T{x) and u^{x). Here we only schematically motivate these results. Intriguingly, 
it turns out that the location of the event horizon rnix^) in the bulk is determined locally 
by the fluid dynamical data at a point (within the derivative expansion), rather than 
globally as usual in general relativity. 

To motivate this physically, within each tube characterized by a given x^, the position 
of the horizon is approximately at r/^ ^ vr T corresponding to that tube. Since T varies as 
a function of x^, so will the horizon position rnix^). In the corrected solution, the surface 
r = 7rT(x) is not the event horizon (for example, it is not a null surface in general), but if 
the variation T{x) is slow, the deviation from the true event horizon is likewise small. 

One can determine the position of the event horizon within our perturbation scheme using 
the fact that the solution settles down at late times to a uniformly boosted planar black hole. 
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In particular, if we expand the horizon location as a series in boundary derivatives, tagged 
as before by e, 

oo 

r = rnix) = tcT{x) + ^ e^r(fc)(x) , 

k=l 

then the coefficient functions are determined algebraically by demanding that the 

surface given by r = rnix) be null. 

A very simple toy model which captures the gist of this argument is given by a time- 
dependent but spherically symmetric black hole, the Vaidya spacetime: 

ds' = - (l - dv^ + 2dvdr + r^ dil^ . 

This metric describes a four-dimensional asymptotically flat black hole accreting null dust, 
so that the mass m{v) increases with time. Assuming that at late times the black hole settles 
down to Schwarzschild, m(f) — m/ as f — )■ oo, and denoting the location of event horizon 
by r = r//(f), we can flnd its position by demanding that it describes the null surface which 
at late times approaches the correct event horizon th — > 2m p as f — )■ oo. Note that the 
normal n to a null surface will be simultaneously tangent to that surface and likewise null. 
For Vaidya, the normal 1-form n = dr — r dv (where ' = ^) is null when 

rniv) = 2 m(f ) + 2 r/^ (t>) rniv) . 

Of course, the exact solution to this equation yields the horizon rniv) non-locally in terms 
of m(f), requiring the knowledge of m(f) for all v < oo. However, when m{v) varies slowly, 
so that m = 0{e), mm = O(e^), etc., we can determine this location in terms of an e 
expansion. For the ansatz 

fuiy) =2m{y) + am{y) m{y) + hm{v) m{vY + cm{vY m{y) + . . . 

this expansion gives a = 8, 6 = 64, c = 32, ... This toy model illustrates that in spite of 
the event horizon being deflned globally (as the boundary of the causal past of the future 
null inflnity) and therefore requiring knowledge of the mass m{y) for all time f < oo, for 
slowly varying m{v) we can nevertheless express rnivo) solely in terms of m its derivatives 
at f = Vq. 

Returning to the problem of interest, we can similarly locate the event horizon rnix^) 
in our dynamical non- uniform planar black hole geometry in terms of T, m'^, and all their 
derivatives, at the given point x^. At flrst order, the position of the horizon is unchanged, 
whereas at second order it is corrected by terms which scale with square of the shear and 
vorticity (see [149], [24] for explicit expressions.) 

Once we identify the position of the event horizon in our geometry, it is easy to check 
that this horizon is regular. In fact, our construction manifestly guarantees regularity: the 
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only curvature singularity of the seed metric is at r = 0, and the source terms which appear 
in correcting the metric order by order do not introduce any additional singularities. The 
final issue to check is the regime of validity of our expansion, and this can be seen to extend 
well inside the event horizon. 

Therefore we have an explicit one-to-one map relating conformal fluid conflgurations 
on R^'^ to asymptotically AdSs inhomogeneous black brane solutions having regular event 
horizons. This is a remarkable statement about gravity, suggesting that fluid dynamical 
conflgurations within the AdS/CFT correspondence naturally uphold Cosmic Censorship. 

To extract further physics from the position of the event horizon, let us consider the 
proper area of its spatial slices. By the second law of black hole mechanics, the horizon area 
cannot decrease with time; or equivalently, the expansion of the horizon generators must be 
non-negative. A well-known identiflcation with thermodynamics translates this statement 
to that of the entropy increasing, or more locally, the entropy current having non-negative 
divergence. Having obtained the event horizon for our geometry explicitly in terms of the 
metric functions u^{x) and T(x), we can verify these statements, and identify the entropy 
current naturally induced on the boundary. 

To obtain the boundary entropy current Jg from the bulk geometry, we can pull-back 
the area form A on the event horizon to the boundary. We perform this pull-back along a 
tube of constant x'^, i.e. along ingoing radial null geodesies. This yields the expression 

+Ci 6/^ + ^2 b'^u^ Vxl" + ... (5. 12) 

with b = ^ and 

1 TT ln2 ,1 ,1 

Ai = - H \ , = — , A = -, 

^ 4 16 4 ' ^ 8 ' ^8' 

Bi = \ , ^2 = ^ , Ci = ^2 = (5.13) 

which can be easily verifled to have non-negative divergence. At leading order, the divergence 
is proportional to a^i, a'^" which is manifestly positive for any non-zero shear. In case of zero 
shear, the relations Bi = 2A^ and Ci + C2 = guarantee non-negativity. We can readily 
see that these relations are indeed satisfled for the entropy current given by our gravity dual 
construction, (5.13). 

The expression (5.12) was written suggestively in a most general form yielding a Weyl- 
covariant entropy current for any set of constant coefficients Aj, Cj. In general, to the 
order we work, the seven coefficients get reduced to five independent ones allowing for Weyl- 
covariant entropy current with non-negative divergence. From the field theory side, this 
would therefore suggest a 5-parameter ambiguity in constructing a sensible entropy current. 
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purely based on the symmetries and the requirement that it correctly reproduces equilibrium 
physics. This is at first sight puzzling, since our gravity construction seemingly fixed all 
these parameters. In fact, this was not quite the case: there is still an ambiguity even on the 
gravity side, corresponding to the freedom to add total derivative terms without changing the 
horizon area, and the pull-back being ambiguous to boundary diffeomorphisms. However, at 
second order this results in a 2-parameter ambiguity for Weyl covariant current with positive 
divergence, so that some mismatch remains. "^^ 



5.4 Generalizations 



The simple analysis presented above was carried out in [22], but this was quickly generalized 
in a number of interesting directions. Here we list some of the most important ones. 



One of the earliest such generalizations involved extending the correspondence to other 
dimensions, relating a ci-dimensional conformal fiuid to asymptotically AdS^+i black 
hole (see [148] for the interesting case of d = 3 and [144, 140] for general^^ d). This was 
already mentioned in §5.3.1 (cf. (5.11)); here we collect the result in more convenient 
form. The most general analysis is presented in [140] who give the second order stress 
tensor (and corresponding transport coefficients), the dual bulk metric along with its 
event horizon and the corresponding entropy current for these fiuid fiows, for a fiuid 
on slowly-varying curved background; see below. In particular, [140] give the stress 
tensor on the d- dimensional boundary metric g^u in a very elegant (and manifestly 
Weyl- covariant) form: 



T, 



flU 



(5.14) 
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dy 



where the pressure p depends on the temperature T and the viscosity of the fiuid rj 
depends on the entropy density s in the usual manner, and t^j denotes a particular 



■^^ Although including higher orders may remedy this discrepancy, [145]. 

■^^ The d = 2 case is rather trivial: the most general solution consists of left and right movers, so the fluids 
cannot locally equilibrate, relatedly the fluid dynamics in 1 + 1 dimensions is equivalent to conservation of 
general conformal stress tensor; also cr^jy = = 0. From the bulk standpoint, any 2 + 1 dimensional 
solution to Einstein's equations with negative cosmological constant is locally AdSs. 
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second-order transport coefficient of the fluid, and 0"^,^ and are the shear and the 
vorticity as before. C^^ap is the Weyl tensor for the metric g^^ in which the fluid lives 
and Vx is the Weyl-covariant derivative. 

One rather intriguing aspect is the striking difference between the phenomenology of 
turbulent flows in 3+1 and 2+1 dimensions, as pointed out in [148]. In the 3+1 di- 
mensional turbulent energy cascade, large scale eddies give rise to smaller scale eddies, 
eventually transferring energy down to scales where viscosity becomes important and 
energy is dissipated. In contrast, the 2+1 dimensional turbulent flows are characterized 
by an "inverse cascade," in which smaller scale eddies merge into large scale eddies, 
creating large long-lived vortical structures. If these qualitative differences extend to 
relativistic fluids, they would suggest a profound difference in gravitational dynamics 
between four and five dimensional gravity. In particular, we might predict that black 
holes in AdS4 would take much longer to equilibrate owing to the fact that the fluctu- 
ations on the horizon could coalesce in macroscopic structures. From the gravitational 
standpoint, this would certainly seem very surprising and counter-intuitive. 

• AdS/CFT, and in particular the fluid/gravity correspondence, extends to bulk space- 
times which are not just asymptotically AdS but also those which are only locally 
asymptotically AdS. In particular, those which asymptote to 

ds^ = ^ [dz^ + dsl^^] (5.16) 

where the 4-dimensional boundary metric ds'^^^ can be any desired slowly- varying met- 
ric. Hence we can consider fluids on curved (and not only conformally flat) manifolds, 
rather than just the Minkowski spacetime R*^"^'^, as has been initiated in [150] and 
carried out to arbitrary dimensions in [140]. It should still remain true that such 
solutions^^ give a universal subsector of CFTs on this curved background, with the 
dynamics given simply by the covariant form of the generalized Navier-Stokes equa- 
tions on dsl^y. As a simple example, [150] (in d = 4) and [140] (for general d) consider 
the explicitly-known case of rotating AdS black hole (which is asymptotically global 
AdSd+i) which describes a rigid fluid flow on S"^"^ x R (see also [151]), con&ming the 
general construction indicated above. 

• In addition, one can include matter in the bulk. This allows for richer dynamics, 
but typically at the expense of losing universality. Early examples of such extensions 
include considering the dilaton (which corresponds to forcing of the fluid) in [150]. 

■^^ For a general (non-slowly varying) boundary metric ds'^^^ we may not know the exact bulk metric 
corresponding to even the global thermal equilibrium. For a discussion when the boundary metric describes 
a black hole with various asymptotics, see [27, 28, 29]; in the high temperature regime the conjectured bulk 
solution, christened 'black funnel', would be a bulk black hole whose horizon stretches all the way out to the 
boundary. 
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• More involved generalizations of the fluid/gravity correspondence with extra matter 
fields include Maxwell U{1) field [152, 153], multiple Maxwell fields and scalars, mag- 
netic and dyonic charges, as well as more exotic models such as those relevant for 
holographic superfluids [154] (see e.g. [24] for further references). 

• Moreover, one can even extend the correspondence to non-conformal fiuids [155, 156] 
as well as to non-relativistic fiuids [157, 158], which allows us to make closer contact 
with familiar everyday systems. 

• Stringy (1/A) and quantum (1/A^) corrections to some of the transport coefficients 
have been computed in [159, 160, 161, 162, 163, 164]. 

Nevertheless, many future directions and puzzles remain, as well as the need for further 
generalizations. For example, of particular current interest is to understand the fiuid/gravity 
correspondence for extremal fiuids which are presently attracting much attention. Also, to 
mimic many of the familiar aspects of fiuid fiows, we need to understand how to confine 
the fiuid within walls in the gravity dual. Still more ambitiously, to understand the rich 
phenomena rooted in quantum processes, we would like to get a better handle on finite- 
effects. 

To summarize, one of the most intriguing features of the fiuid/gravity correspondence is 
that it provides us with a window into the generic behavior of gravity in a nonlinear regime, 
mapping long-wavelength (but arbitrary amplitude) perturbations of AdS black holes to the 
more familiar physics of fiuid dynamics. Apart from the obvious conceptual advantages, 
one has a tremendous computational simplification for numerical studies of gravitational 
solutions since the fiuid dynamics lives in one lower dimension. 

The bulk spacetime solutions discussed here describe a generic (time-dependent and non- 
uniform) planar black hole. Each of these solutions (with regular fiuid data T[x) and Ufj_{x)) 
has its singularities hidden from the boundary by a regular event horizon. In this sense, all 
gravitational solutions dual to regular solutions of fiuid dynamics uphold the cosmic cen- 
sorship conjecture. At the technical level, the key to the fiuid/gravity construction lies in 
utilizing the long-wavelength regime of fiuid dynamics. This allows effectively a linearization 
of Einstein's equations in that one ends up solving linear ordinary differential equations to 
solve for the bulk metric, and corresponding boundary stress tensor, to any order in a bound- 
ary derivative expansion in a completely procedural manner. Furthermore, having obtained 
the bulk geometry, this expansion likewise enables us to determine the radial position of 
its event horizon, remarkably expressed locally in terms of the fluid dynamical data. One 
can then easily verify that the area of this event horizon is necessarily non-decreasing, in 
accordance with the Area Theorem. The corresponding entropy current induced from the 
bulk solution is then guaranteed to satisfy the Second Law of thermodynamics. 

The boundary fluid stress tensor contains new quantities of interest, namely the various 
transport coefficients which characterize the fluid. At flrst order one recovers the previously- 
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known value of viscosity and verifies that ^ = More importantly, the construction allows 
the extraction of second order fluid parameters (r^r, Ai, A2, A3) (see also [83]). This has been of 
interest in QCD phenomenology, especially in understanding certain characteristic features 
of the quark-gluon plasma. 

6 Beyond long- wavelength approximation 

While impressively powerful, the fluid/gravity correspondence as formulated is only valid 
in the long-wavelength regime. Albeit useful to study gravitational duals of fluid flows, it 
falls short of describing more interesting aspects of dynamics, viz., the approach to local 
equilibrium and the behaviour of general fluctuations which are not necessarily of long- 
wavelength. Likewise, from the bulk point of view, the long-wavelength regime does pose 
a substantial restriction, in the sense that many 'natural' bulk solutions do not fall in this 
regime. 

The simplest example is the static spherically symmetric Schwarzschild-AdS black hole 
in global AdS, whose radius is larger, but not parametrically larger, than the AdS radius, 
r+ > -RacIS- Such large (ish) global Schwarzschild-AdS black hole does not fall within the 
long-wavelength regime because the boundary curvature is comparable to the thermal scale 
(for instance the term proportional to the Weyl tensor in (5.14) obviates the expansion). 
However, this is a regular solution, admissible from the point of field theory, describing 
simply an equilibrated plasma on the boundary sphere. Furthermore, one expects (assuming 
cosmic censorship) that this solution actually captures the end-point of a generic collapse 
scenario in AdS for initial data whose energies are in the appropriate range. This then 
suggests that the dual state in the CFT is not described by a fluid stress tensor in spite of 
being globally equilibrated and exhibiting no spatial variations! 

One would therefore like to understand the generic situation where perturbations take 
us away from hydrodynamic equilibrium, involving excitations of high frequency modes. A 
related question of interest concerns the evolution of a highly excited state in a quantum 
system. Consider preparing a quantum state, with a large vacuum expectation value for 
the Hamiltonian operator, which is not an energy eigenstate, as can be done for instance 
by injecting energy into the system. One expects this state to explore the available phase 
space and evolve dynamically to a thermal density matrix, in the sense that the late time 
observables are thermal expectation values. In such situations one would like to track the 
evolution of such a state. From the bulk standpoint this evolution should be associated 
with the formation of a black hole, and questions like the approach to equilibrium and 
validity of hydrodynamics can be addressed once the solution is understood. To summarize 

""^This concept is also of interest in condensed matter systems, where it is referred to as quantum quench. 
See [165, 166, 167, 168] and references therein for details on this subject. We should also note that there has 
recently been some progress in holographic modelling of quantum quench behaviour in [169]. 
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our motivation, we would like to move away from the axes in Fig. 1 and explore the region 
of arbitrary amplitudes and momenta. 

In this section we will review various works that have addressed these issues. While most 
of the techniques for investigating the general evolution involve numerical simulation (we are 
in the non-linear regime with no control parameter generically) , there are some interesting 
examples where analytic techniques can be brought to bear. We will begin by describing 
simple toy models which allow detailed exploration of time dependence and then discuss 
some concrete attempts to study black hole formation in the AdS/CFT context. Towards 
the end we revert back to the story of equilibrium dynamics, but for field theories in curved 
spacetime. 

6.1 CFTs with strong time dependence 

We have seen above that even innocuous static configurations can lie outside the long- 
wavelength regime, by virtue of living on a compact space. However, of greater interest in 
the present setting are configurations exhibiting explicit time-dependence, which are typically 
outside the long-wavelength regime. This is because a typical dynamical process involving a 
black hole (such as its formation or coalescence, quasinormal ringing, etc.) occurs on time 
scales set by the black hole size, which in the regime of interest is given by the thermal scale, 
rather than being parametrically longer. Nevertheless, there are interesting and important 
physical processes which may be at least partially understood from the fluid/gravity frame- 
work. In the following discussion, we will first consider systems with 'mock' time-dependence 
given by simply considering a static configuration in different coordinates, and then discuss 
several configurations which are genuinely time-evolving. We will see that even the first, 
seemingly rather trivial, class offers intriguing surprises. 

6.1.1 Mock time-dependence with evolving event horizon 

One of the simplest examples of a secretly static configuration which nevertheless provides 
an instructive toy model for a time-dependent scenario is given by the conformal soliton 
geometry originally introduced in [132], and then studied extensively in [170]. The spacetime 
in question is simply a patch of the well-known Schwarzschild-AdS black hole; so the explicit 
metric is known exactly, and admits a high degree of symmetry. Nevertheless, if we work in a 
coordinate system such that the field theory lives on fiat space R^'^ rather than the Einstein 
Static Universe x R^, which corresponds to considering only a 'Poincare patch' of the 
Schwarzschild-AdS black hole, the time translation symmetry is no longer manifest, and 
the solution looks highly dynamical. Pictorially, it describes a black hole entering through 
the past Poincare horizon and exiting through the future Poincare horizon, with its closest 
approach to the boundary occurring at t = (Poincare time). In the boundary CFT, this 
corresponds to a finite energy lump which collapses and re-expands in a time-symmetric 
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Fig. 5: Conformal Soliton energy density (left) and dual geometry (right). On left, TQQ{t,r) is plotted 
as a function of boundary time t and boundary radial variable r. On right, the outer cylinder 
represents the AdS boundary, the inner (red) cylinder corresponds to the global Schwarzschild- 
AdS event horizon, which coincides with the apparent horizon of the conformal soliton, whereas 
the flared (blue) surface is the actual event horizon. 

fashion (cf. the left panel of Fig. 5). Here the fluid dynamical approximation is not valid at 
all times, but because this fluid flow is conformal to a stationary fluid on the Einstein static 
universe, the stress tensor is shear-free; that is, there is no dissipation in this fluid flow. 

Nevertheless, the event horizon in this case is rather interesting. Since our coordinate 
patch includes only part of the future infinity of the global Schwarzschild-AdS, the actual 
event horizon for the conformal soliton lies outside the global event horizon of Schwarzschild- 
AdS. For orientation, the event horizon for both global Schwarzschild-AdS and the conformal 
soliton is plotted in right panel of Fig. 5. In fact, [170] observed that the area of this conformal 
soliton event horizon diverges at late times. This surprising result led to a puzzle: if we 
associate the entropy of the corresponding CFT conformal soliton state to the area of the 
conformal soliton event horizon, as is usually assumed to be the case, then we find that 
this entropy likewise diverges at late times. But the conformal soliton describes a shear-free 
flow, with no entropy production whatsoever. Said differently, it is easy to see [170] that the 
conformal transformation from the CFT on x to the CFT on R^'^ leaves the entropy 
invariant. But the former describes a perfect fluid in global thermodynamic equilibrium: its 
entropy is finite and constant in time. 

Already prior to [170], it has been argued in several different contexts [171, 172] that it 
may be more appropriate to associate the entropy of the CFT configuration to the area of 
the apparent horizon, rather than the event horizon, in the bulk dual. In fact, as explained 
in [170], the apparent horizon in Poincare slicing coincides with the global Schwarzschild- 
AdS event horizon, whose area is indeed constant and given by the expected value. Since the 
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event horizon and the apparent horizon behave radically differently (cf. Fig. 5), this conformal 
soliton geometry provides a good testing ground for studying the distinction between the 
two horizons and the role they play for the associated CFT dual. We see that in this case 
the CFT entropy is clearly more naturally associated with the apparent horizon rather than 
the event horizon. 

The large disparity between locations of the event and apparent horizons might seem 
enigmatic in light of [149], where it was argued that within the fluid dynamic regime, the 
event horizon and the apparent horizon track each other closely. The essential difference 
between these two contexts is that in [149], it was assumed that the geometry would settle 
down at late times to a stationary finite-temperature black hole. This is not the case for the 
conformal soliton geometry. It is precisely these late time boundary conditions that force 
the apparent horizon and the event horizon to behave very differently. 

This observation might prompt one to expect that if the spacetime does not settle down 
at late times, then the entropy will not be given by the event horizon area. However, we will 
now briefly mention another example, this time with genuine time dependence, where this is 
not the case. In particular, let us consider the holographic dual to Bjorken flow [173, 174], 
also studied in [170]. Since Bjorken flow plays a central role in understanding the post- 
thermalization evolution of the QGP produced in heavy-ion collisions, this is a more broadly 
useful example to study. The basic physical picture developed in [175] is that, assuming local 
thermal equilibrium, in the central rapidity region of ultra-relativistic collisions of heavy ions 
one can model the flow of the plasma via quasi-ideal hydrodynamics. In the hydrodynamic 
description, it is assumed that the fluid evolution respects the boost symmetry along the 
collision axis. This implies a boost invariant expansion of the fluid, consistent with the 
observed distribution of the particles in the collision process. 

In [173, 174], the corresponding bulk geometry was constructed as a perturbation expan- 
sion in the boundary time coordinate which is valid at late times. By demanding regularity 
of the solution at leading orders, the authors were able to derive the transport properties of 
the plasma, in particular the shear viscosity r] which saturates the bound r]/s > I/Att [15] 
mentioned in §4.2. The study of the gravitational dual at higher orders was undertaken to 
derive the relaxation time of the plasma in [176, 83]. Although initially the regularity of 
the dual spacetime was brought into question as subleading singularities were encountered 
[176, 177], this issue was subsequently addressed in [178, 179] where the authors used the 
framework of the fluid/gravity correspondence [22] to argue that the spacetime was indeed 
regular. Explicit confirmation was provided in [170] by constructing the global event horizon 
for these geometries, finding that it is regular and closely tracks the location of the apparent 
horizon, as previously computed in [179] and verified in [170]. 

We will revisit the intriguing question of whether this can be taken as a general tenet in §8. 
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6.1.2 Black hole formation in AdS and approach to equilibrium 

Having discussed in §6.1.1 a simple case where one engineers apparent time dependence by 
a non-trivial coordinate change in a static solution, we now turn to recent discussions of 
explicit time dependent phenomena in AdS/CFT. These examples provide a clean way to 
understand the deviations from equilibrium and approach to the same from an arbitrary 
initial state. 

One way to study time-dependent phenomena in AdS/CFT is to start with a system 
in the vacuum state and inject in energy at some suitably chosen moment, say t = 0. At 
early times t < the system is holographically described by the pure AdS spacetime, while 
after the perturbation one needs to solve the bulk Einstein's equations in order to infer the 
nature of the evolving geometry. The injection of energy can be achieved by perturbing the 
boundary Hamiltonian with an explicit insertion of a source at t = 0. In the bulk description, 
the source corresponds to a change of boundary conditions for the corresponding field; in 
effect one has to turn on non-normalizable modes for the bulk fields. Equivalently, one can 
consider sending in a null shell at t = from the boundary, which deforms the spacetime 
from pure AdS in the interior. In either case we expect to seed a black hole collapse; the 
energy contained in the deformation (or in the null shell) will cause a black hole to form 
in the interior of the spacetime, which at late times will asymptote to a new steady state 
solution. A schematic Penrose diagram of such a process is sketched in Fig. 8. 

In [172] the first non-trivial such calculation was carried out.^'^ The idea was to inject 
energy into the boundary field theory by changing the boundary metric. Since the boundary 
metric couples to the energy-momentum tensor, we are effectively pumping energy into the 
system directly. The set-up used by [172] was to consider a background for the field theory 
which asymptotes to the fiat Minkowski background in the far past and far future, but with 
different length scales. In particular, they consider explicitly time dependent background of 
the form: 

rfsLy = -dt' + e^»W rfxl + e-2^°(*) dxl . (6.1) 

With a suitable choice of Bo{t), e.g. approximating a step function obeying the early and 
late time boundary conditions, one can inject energy into the system as a consequence of 
the time- varying spatial length scales. Physically the change in length scales corresponds to 
anisotropizing the pressure in the various directions. 

Since the deformation is purely in the metric degrees of freedom, the problem can be 
treated within the consistent truncation of pure gravity in AdS^+i, i.e., by solving (2.6) for 

One can of course construct such geometries by patching together spacetimes across thin domain waUs 
using the Israel junction conditions, but since we are interested in real dynamics we need to actually solve 
the bulk equations of motion. 

'*'^In fact earlier analysis of equilibration in strongly coupled systems was considered in [180] building on 
the construction of time dependent Janus type geoemtries in [181]. 
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the bulk metric Gmn subject to the boundary condition of approaching (6.1) asymptotically. 
The numerical solution of this problem clearly illustrates the process of black hole formation; 
in particular, [172] were able to demonstrate the presence of an event horizon by examining 
the causal past of the AdS boundary. Moreover, the computation of the boundary energy 
momentum tensor using the holographic recipe (2.16) showed clearly that the pressure in the 
system starts out being anisotropic between the longitudinal (xy) and transverse (xj_) direc- 
tions, but evolves towards an isotropic equation of state characterizing eventual approach to 
equilibrium. In [182], the phenomenologically interesting boost invariant Bjorken flow [175] 
was also modeled in the bulk by a numerical solution. This calculation provides a useful 
way to make contact with the hydrodynamic description of the system at late time for which 
approximate solutions in the fluid/gravity sense were constructed in [173, 174] as described 
at the end of §6.1.1."^^ Numerical investigations of the early time behaviour relevant for the 
Bjorken flow have also independently been considered in [184]. 

While the constructions described above involve numerically integrating Einstein's equa- 
tions, it is nevertheless possible to understand aspects of black hole formation and ther- 
malization in AdS spacetimes analytically. This was achieved in [185] by considering a 
perturbative solution in the amplitude of the fluctuation. The basic idea is to start with the 
initial data similar to the one described above, but to take the injection of energy to be of 
small amplitude to allow for a perturbation expansion in amplitude. The key ingredient 
in the analytic construction was the observation that the perturbation series in amplitude 
needs to be resummed, for the naive expansion fails to converge. This is perhaps not so sur- 
prising given that a black hole, no matter how small, is necessarily a large perturbation on 
causally-trivial constant curvature spacetime. The resummed perturbation series naturally 
leads to a change in the background reflecting the process of black hole formation. Assuming 
that one started with pure AdS spacetime and injected energy (or sent in a null scalar pulse) 
at t = 0, one flnds that the resummation of the perturbation series shifts the background 
towards a Vaidya type solution, which describes gravitational collapse. 

The analytic results of [185] illustrate two key features of gravitational collapse in AdS 
spacetimes: 

• In global AdSd+i spacetimes there is a threshold for black hole formation if one exam- 
ines the system on a time-scale of the order of the AdS scale after the perturbation. 
This is similar to the gravitational collapse in asymptotically flat spacetimes where 
the space of initial data splits into distinct domains separated by a co-dimension one 
critical surface which captures the threshold conflgurations for black hole collapse. 

• Arbitrarily small energy densities lead to black hole formation when the boundary is a 

"^^ These numerical techniques have also been adapted to study the evolution of perturbations on unstable 
phases of the field theory, see [183]. 

'^^ In addition to working in pure gravity, [185] also considered the somewhat simpler context of a scalar 
collapse in AdS. 
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(deformed) planar BH 
(a) (b) 

Fig. 6: A sketch of our two novel classes of solutions dual to a deconfined phase of a strongly coupled 
field theory on a black hole background: (a): black funnel and (b): black droplet above a 
deformed planar black hole. The shaded part indicates regions which are behind horizons in 
the bulk spacetime, and the thick (red) line indicates the black hole of the boundary metric. 



copy of Minkowski space; in this context there is no threshold for black hole formation. 
6.2 CFT on non-trivially curved backgrounds 

In §6.1.2 we indicated how one might perturb a system in equilibrium by deforming the 
spacetime background in which it lives. While this is interesting in order to understand 
deviations from equilibrium, examining the behaviour of quantum fields in curved spacetime 
brings another dimension to the problem. As described towards the end of §2.1 the AdS/CFT 
correspondence provides a tool to investigate such issues as well. In general one needs to 
construct all bulk manifolds whose boundary is the background on which the field theory 
dynamics takes place. 

Moreover, from the standpoint of our general discussion, static configurations might 
have explicit non-long-wavelength feature resulting from the boundary background geometry 
the CFT lives on. For example, consider the field theory on a Schwarzschild black hole 
background. Far away, the equilibrium state is simply a thermal state at a temperature 
given by the black hole's temperature, which in turn determines the microscopic scale £mfp ~ 
l/^BH- But the curvature of the background geometry is of that same scale near the black 
hole. From the dynamical point of view, the fluid near the black hole wants to fall in on a 
timescale determined by the black hole, i.e. also ^mfp- Considering such system in equilibrium 
with the outgoing Hawking quanta leads to an interesting bulk dual: [27] suggested that the 
dual geometry exhibits new (yet to be found) solutions, dubbed black funnels and black 
droplets. For ease of orientation, these are sketched in Fig. 6. Far away from the black 
hole in the boundary directions, the bulk configuration is well-described by the fluid/gravity 
framework, approaching the planar black hole at the requisite temperature (namely the local 

As a follow-up, evidence for lower-dimensional realization of these was analyzed in [28] and further 
discussion of field theory on AdS black hole background was presented in [29] . 
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temperature of the boundary field theory heated up by the boundary black hole). However 
the event horizon of the background spacetime which the field theory lives on must be 
continued into the bulk on the bulk event horizon. Near the boundary, such an extension 
looks like the AdS black string, but in the bulk this black string is deformed. If it joins with 
the bulk planar black hole and forms a single connected horizon, the solution is called a black 
funnel; on the other hand, if it caps off and only slightly deforms the disconnected planar 
black hole horizon, it is called a black droplet. Depending on which geometry is realized, 
the physics (such as the response of the system to perturbations) differs significantly [27]. 
As the parameters of the background are varied, we furthermore expect a phase transition 
to occur. Such rich physics would of course not be apparent had one attempted to use the 
fluid/gravity framework with thermal fluid stress tensor everywhere (for example we could 
never see the droplets this way), but then one would be alerted to being well out of its regime 
of validity by the description breaking down. 

These examples serve to illustrate the rich physics exhibited by systems which are not 
in local thermodynamic equilibrium but which are nevertheless accessible for exploration 
within AdS/CFT. 

7 Diagnostics of dynamics 

Let us pause to take stock of the discussion so far. In §3 we have discussed how to extract 
certain properties, such as transport coefficients of the CFT by considering infinitesimal 
deviations from equilibrium, which is fortuitously captured by probes in the equilibrium 
configuration itself. We then went on to consider genuinely time-dependent configurations 
in §5, and explicitly constructed the bulk spacetime and the corresponding CFT fluid stress 
tensor to second order. While the deviations from global equilibrium can be large, the spatial 
and temporal dependence had to be sufficiently slow for the boundary derivative expansion 
to be valid. This posed a rather severe restriction, as many processes of interest happen on 
faster time scales, comparable with the thermal scale. Hence in §6 we tried to venture beyond 
this long-wavelength regime. Except for toy models with 'mock' time-dependence, we have 
much less handle on such systems, but these strongly time-dependent cases nevertheless yield 
more insight about generic physical processes. 

In such genuinely time-dependent scenarios, it is often the case that we have a good 
idea of what the bulk spacetime looks like - typically we know the explicit metric, at least 
approximately, at least in some region. However, having handle on the bulk does not mean 
that we know much about the CFT dual; actually constructing the CFT state dual to the 
given bulk configuration is well beyond our means at any rate, though we can use the 
gravity side calculation to read off various observables in that state. For example, using the 
asymptotic fall-off of the metric, we can determine the stress tensor expectation value in the 
dual CFT state. Conversely, we can turn the logic around and ask: which CFT data should 
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we use to learn as much as possible about the bulk spacetime? 

The purpose of this section is to review some of the CFT 'observables' which reveal the 
salient features of the bulk geometry, focusing on probes which we can use in time-dependent 
cases. The philosophy here is to suppose that we have at our disposal all that we wish to 
know about the CFT state, and ask where should we look to see a particular bulk feature. 
As already indicated in the Introduction and reviewed in §2, this is not an easy question, 
since bulk locality is not manifest, nor indeed well-understood, in the dual CFT. On the 
other hand, this was part of the motivation of the works reviewed below: by determining 
the signature of specific local bulk features in the CFT dual, we gain some insight into how 
holography encodes bulk locality, and in turn into the emergence of spacetime (see [186, 187] 
for recent progress). 

In classical general relativity, concepts such as causal structure, event horizons, sin- 
gularities, etc., play an important role in understanding the spacetime geometry. In the 
semi-classical approximation these concepts are useful in understanding the dynamics of 
quantum fields in curved backgrounds. Given the AdS/CFT correspondence, it is therefore 
interesting to understand the field theoretic encoding of these geometric features. From the 
central role played by geometry in classical general relativity, one naively expects them to 
have a well defined representation in the field theory. For example, consider a time depen- 
dent process of black hole collapse in the bulk. Since the formation of the event horizon is a 
sharply-localized spacetime event (its teleological nature notwithstanding), we would expect 
that this will manifest itself in some correspondingly sharp feature in the gauge theory. This 
expectation was indeed verified in [188], as we explain below. Of course, if we could extract 
the full spacetime metric from the CFT data, we could then reconstruct these salient fea- 
tures of the geometry. However, that is not a very direct route, since extracting the actual 
spacetime metric is harder than extracting its geometrically interesting features. 

One might naively think, using the ideas of §5, that given the stress tensor expectation 
value everywhere, one can nevertheless reconstruct the entire bulk metric. This is not cor- 
rect. The stress tensor expectation value only knows about the asymptotic behaviour of the 
metric. In case of holographic RG [26], one can write the metric in a radial series expan- 
sion around the boundary, but there is no guarantee that this series will converge inside the 
bulk. Indeed as explained in [22], a general (non-fluid) conformal stress tensor will lead to 
naked singularities in the bulk. On the other hand, within the fluid/gravity correspondence 
framework reviewed in §5, we can reconstruct the spacetime exactly in the radial direction 
well inside the (regular) event horizon - but only at the expense of confining ourselves to 
the long-wavelength regime. In other words, in such a class of spacetimes, the asymptotic 
behaviour determines the geometry in the bulk as well, in the 'tubewise' manner described 
above. While it is impressive that local expectation values carry so much information about 
the bulk^"^, it is nevertheless clear that we will need at least bi- local observables, such as 

This is intimately tied to the AdS asymptotics, which in a sense 'magnify' the signals to be discernible 
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correlation functions, to learn about the full bulk geometry. 



7.1 Extracting bulk metric from CFT correlators 

The conventional lore concerning the UV/IR relation, which maps local regions in the in- 
terior of the bulk to non-local objects in the CFT, naively appears to preclude extracting 
much useful information about bulk geometry, especially any precise signal about the causal 
structure in the neighbourhood of the event horizon. The initial clue on how the field theo- 
retic observables transcend the classical barrier of the event horizon to encode information of 
behind-the-horizon physics^^ emerged from the progress made in identifying the CFT signa- 
ture of the black hole singularity [191, 192], which built on the prior work of [193, 194, 195]. 
Schematically, the idea was to use the intrinsic non-locality of the boundary correlation 
functions to identify signals of the bulk curvature singularity. 

The basic strategy was to look at the bulk Green's functions in a saddle approximation 
where they are dominated by geodesies and thereby use simple bulk computations to extract 
the behaviour of the boundary correlation functions. When spacelike geodesies which get 
repelled by the curvature singularity become almost null, they provide a large contribution 
to the correlation function, akin to light-cone singularities, in the appropriate regime. At 
the same time, such geodesies penetrate arbitrarily close to the curvature singularity, where 
their contribution is most significant. Hence the black hole singularity is manifested by 
a particular light-cone-like singularity in the field theory correlation function. The left 
panel of Fig. 7 sketches the relevant spacelike geodesies for the original Schwarzschild-AdSs 
spacetime; as a critical time tc is approached by the geodesic endpoints (which we can 
think of as the correlator insertion points), the geodesic approaches two null geodesies which 
meet at the singularity. Moreover, it has been proposed in [196] that this technique can 
be exploited to extract physics from further beyond the event horizon, and in particular 
used to investigate aspects of inflationary geometry within AdS/CFT. On the right panel, 
such almost-null geodesies are indicated for the far more complicated spacetime considered 
in [196] containing a de Sitter universe behind the horizon of a Schwarzschild-AdS black 
hole. Intriguingly, unlike for the previous case of just the Schwarzschild-AdS geometry, in 

from the boundary. A very early example which illustrates that far more information about the bulk is 
contained in the local CFT expectation values than naively expected from the UV/IR correspondence is 
given in [189], where the authors show that even such detailed information as the size of a small i?Ads) 
object can be read-off from local CFT expectation values. 

*® One of the early motivations that this should be possible in the first place, despite bulk causal obstruc- 
tions, was presented in the gedanken-experiment of [190]. 

The identification of the singularity in the strict large N limit is easiest when formulated in terms of 
momentum space correlators [192] as opposed to direct computation in position space where the singularity 
is not visible in the primary sheet of the correlator [191]. In particular, it was confirmed in [192] that the 
signatures of the singularities disappear for finite rank of the boundary theory gauge group, implying that 
the singularities are resolved in quantum gravity. 
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Fig. 7: Illustration of mapping out bulk causal structure by analytically continued correlators, as 
discussed in [191] for the Schwarzschild-AdS5 spacetime (left panel) and in [196] for a more 
complicated geometry describing deSitter universe behind the horizon of a Schwarzschild-AdS 
black hole (right panel). In both Penrose diagrams, the vertical solid lines represent the AdS 
boundaries, the top and bottom (purple) curves the curvature singularities, the dashed diagonal 
lines the horizons; on right the horizontal lines are de Sitter boundaries, the vertical dashed 
line the de Sitter origin, and the vertical (brown) curve a domain wall separating the dS and 
SAdS regions. All other curves correspond to the geodesies discussed in the text. 

these more complicated cases tlie relevant correlator insertion points appear on the same 
boundary. 

As evident, these correlators are therefore highly sensitive to certain bulk information 
which might otherwise be completely inaccessible to a local bulk observer. In fact, even such 
subtle differences as the fuzzball picture of the black hole [197, 198] and a genuine eternal 
black hole geometry would produce glaringly different signatures in these correlators [199]. 
For instance, as discussed in [200] the fact that the eternal black hole corresponds to a den- 
sity matrix can be used to argue that the correlation functions computed from this geometry 
should be periodic in imaginary time. This is of course well known from the discussions of 
computing thermal correlation functions [56, 57, 62]. However, any single horizon-free geom- 
etry will give rise to correlation functions that do not display this periodicity in imaginary 
time. 

This geometric picture can therefore be used to distinguish black hole spacetimes from 
the putative horizon-free microstate geometries. To understand this issue further recall that 
while the computation of the Green's functions in the black hole geometry leads to field 
theory observables in the canonical (or more generally grand canonical) ensemble, the com- 
putation in a single microstate geometry leads to a pure state correlator. One can of course 
consider an ensemble of pure states and indeed one would then recover the thermal answer 
via appropriate ensemble averaging - this then corresponds to a micro-canonical ensemble 
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computation, which by ensemble equivalence should agree with the thermal result, as was 
convincingly argued in [201, 202]. However, as demonstrated in [199], by examining the 
detailed analytic structure of the correlation functions, one can amplify the small differences 
between the micro-canonical and canonical computations and use this as a distinction to tell 
apart spacetimes with and without horizons. 

Given the remarkable level of detailed bulk information encoded in such CFT correlators, 
one might wonder how far-reaching this tool is for 'decoding the hologram' in general situa- 
tions. The main shortcoming of this method is that it is most useful for static^° geometries, 
where we can determine the correct physics using the Euclidean continuation. Nevertheless, 
this limitation can be overcome, as demonstrated by [188], where the authors consider a 
manifestly time-dependent spacetime describing gravitational collapse and argue that the 
horizon formation event can be detected in the boundary theory by examining the structure 
of singularities of generic Lorentzian correlators. The basic idea is that CFT correlators 
will exhibit light-cone singularities when the operator insertion points are connected by a 
strictly null geodesic. Such a geodesic may be confined to the boundary (in which case 
the corresponding light-cone singularity is a familiar feature of the field theory), but more 
interestingly it may also penetrate into the bulk spacetime. In the latter case, the connec- 
tion implies that CFT correlators in excited states have additional Lorentzian singularities 
inside'^ the boundary light cone, which in [188] were christened hulk-cone singularities. 

This remarkable result has important and useful applications. Because the endpoints of 
null geodesies through the bulk depend on the bulk geometry, the locus of the bulk-cone 
singularities changes as the bulk geometry changes. In [188] the nature of these changes 
was explored for AdS deformed by a presence of a radiation "star", a collapsing shell, and 
an eternal black hole; in each case, the pattern of bulk-cone singularities exhibits quali- 
tatively distinct features. In particular, in the collapse scenario, [188] demonstrate that a 
sharp horizon-formation time can indeed be extracted from the pattern of singularities; cf. 
Fig. 8: Consistently with the teleological nature of the event horizon, the relevant bulk cone 
singularity occurs for characteristic insertion points, one at infinite time t = oo, the other 
null-separated from the location of the horizon formation event, t = th- This is because the 
corresponding null geodesic (drawn as diagonal blue dotted/dashed line from th in Fig. 8) is 
the latest null geodesic to reach the boundary. Note that the time th precedes the formation 
of the shell, tg, while tn and ts are spacelike-separated. Moreover, as the black hole is about 

Following conventional terminology, by static we mean admitting a Killing field which is asymptotically 
timelike, but not necessarily globally static. In particular, this includes nontrivial geometries such as black 
holes or de Sitter, which have regions of strong time-dependence. Indeed, it is precisely due to these regions 
that the CFT signals proposed in [191, 196] are present. 

As proved by [203], for certain wide class of spacetimes (including the present case of interest) with 
timelike conformal boundary, any fastest null geodesic connecting two points on the boundary must lie 
entirely within the boundary. The vacuum state of the CFT plays a distinguished role, in that all null 
geodesies through pure AdS bulk take equal time to reach the antipodal point as the boundary null geodesies. 
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Fig. 8: Penrose diagram for collapse in AdS. The vertical solid line represents the AdS boundary, 
the vertical dashed line the AdS origin, and the upper wavy curve the curvature singularity, 
produced by a null shell which originates on the boundary at time ts (CFT dual of this process 
is discussed in §6.1.2), and collapses to zero area. In the bulk, event horizon (dashed diagonal 
blue line) nucleates at bulk time tjj, whose signal can be seen by bulk-cone singularity with 
one insertion at and the other at t = oo. 



to form, the separation (to — ti) between the insertion points (ti, to) of the bulk-cone singular- 
ities grows without bound as U ^ ; in fact, the exact evolution of this separation should 
carry information about criticality in black hole collapse, such as that present in context of 
Choptuik scaling [204]. 

Given that distinct geometries typically lead to distinct endpoints of null geodesies 
through these geometries, one might well ask whether one can invert this relation to ac- 
tually extract the bulk metric from the locus of the bulk-cone singularities. At first sight this 
appears a daunting task since each geodesic passes through a 1-dimensional set of spacetime 
events, so there is no natural map between a single geodesic and a single spacetime event. 
On the other hand, each spacetime event in {d + l)-dimensional bulk is pierced by a (c? — 1)- 
dimensional set of null geodesies, each specified by its direction at that event. Moreover, by 
continuity nearby geodesies pass through nearby events; but since the separation between the 
nearby geodesies is typically larger in some regions than in others, such regions dominate the 
corresponding changes in the bulk-cone singularities. This motivates the expectation that 
if one takes the combined information carried by an appropriate set of geodesies, one can 
recover information about the metric in a localized region. However, to what extent one can 
carry out this metric extraction in general remains an interesting open question. 

The simplest setting in which to investigate this question is to focus on static, spherically 
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symmetric spacetimes. The most general such spacetime can be written in the form 



where /(r) and h{r) are arbitrary^^ functions of one variable with the prescribed boundary 
conditions. Spherical symmetry allows us to consider only the equatorial plane, so that the 
geodesies can be described in 3-dimensional (t, r, ip) space, and from each boundary point 
[t, 99) at r = cxD there emanates a 1-parameter family of null geodesies, specified by the 'an- 
gular momentum per energy' i. On the boundary, time translation and rotational invariance 
imply that the bulk-cone singularity locus is specified by {At{i), Aip{i)), i.e. one function's 
worth of data At{Aip). It is then natural to expect that without additional information, this 
single function does not suffice to recover two independent functions /(r) and h{r). Indeed, 
from the bulk perspective, this is consistent with the fact that null geodesies are insensitive 
to conformal rescaling of the metric. Nevertheless, the above counting suggests that once 
we fix the conformal factor, we should be able to recover the metric in at least some region. 
Explicit reconstruction was carried out in [205, 206] for spacetimes of the form (7.1) with 
h{r) = l//(r), using numerical and analytical means respectively. Considering various ad- 
hoc examples of f{r), the authors demonstrated that using only the corresponding bulk-cone 
singularities, /(r) can be recovered for all r > r^i^, where either r^^^ = or it corresponds to 
a null circular orbit. Furthermore, if the available boundary data includes the endpoints of 
not only null, but also spacelike geodesies at constant t, these works demonstrated extraction 
of both /(r) and h{r) in (7.1). 

These examples illustrate that by using the non-local nature of correlation functions one 
can obtain information about the details of the bulk spacetime. In fact, this is already a vast 
improvement to the previously-mentioned avenues to decoding the bulk; in particular, it is 
not reliant on long-wavelength regime. Moreover, in contrast to the previous discussions of 
probing the black hole singularity in AdS/CFT [191, 192], the use of bulk-cone singularities 
does not rely on analytic properties of the CFT correlators. Nevertheless, as indicated above, 
bulk cone singularities have several shortcomings of their own, primarily related to general 
properties of null geodesies. First, they are insensitive to conformal rescaling of the bulk 
geometry; hence they cannot determine the metric fully. Second, since null geodesies are 
causal, bulk cone singularities can never be used to probe the geometry inside an event 
horizon. In fact, using the methods mentioned above for recovering /(r) = l/h{r) in (7.1), 
one cannot even probe the geometry beyond a null circular orbit (which is separated from the 
horizon by a factor of (9(r+)). Finally, another limitation of using the bulk-cone singularities 
to extract information about the bulk geometry arises from the fact that it is well-tailored 
mainly to asymptotically globally AdS spacetime, i.e. the field theory living on a compact 
space. When the field theory lives on e.g. R^'^, there may be no bulk-cone singularities, as 

Here we ignore the bulk field equation, or equivalently assume no knowledge of the bulk matter content, 
and ask whether the metric can be recovered solely from the restricted set of boundary data given by the 
bulk-cone singularities. 
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demonstrated in [188] for any state respecting the same symmetries. 



7.2 Detour: surfaces of different dimensionality as probes of geometry 

So far, we have seen that although CFT correlators, in particular the locus of bulk-cone 
singularities, may encode a large amount of information about the bulk geometry and for 
certain class of states recover the bulk metric entirely, they nevertheless do not suffice in 
general to recover the full metric. We will therefore take a brief detour to motivate our 
expectations as to which CFT observables would be likely to provide better probes of the 
bulk geometry. 

Typically, geometrically defined objects such as geodesies, extremal surfaces, etc., which 
are anchored on the boundary, correspond to some probe (or observable) in the dual CFT. 
For example, we have already seen the relation between geodesies and correlators in §7.1, we 
have mentioned the well-known relation between the area of two-dimensional minimal surface 
describing a string world-sheet and the Wilson-Maldecena loop [124] in §4.3.1, and in §7.3 
below we will use the relation between co-dimension two extremal surface and entanglement 
entropy [207, 208]. One may therefore ask, purely as a geometrical question in the bulk, how 
far into the bulk can such surfaces penetrate. The motivation behind this question is the 
tacit assumption that we can extract the information about the geometry to wherever the 
bulk dual of our CFT probes can reach. This is in turn motivated by our experience with 
geodesies mentioned above, as well as with certain higher-dimensional surfaces discussed in 
§7.3 below. Conversely, it is clear that if the bulk dual of a given CFT probe does not pass 
through a specific bulk region, then it does not carry any information about the spacetime 
geometry in that region. 

Let us for simplicity continue considering only static spherically symmetric spacetimes of 
the form (7.1) discussed above. Geodesies in such a spacetime can be described in terms of 
effective potential describing the radial motion and conserved quantities E and L conjugate 
to the temporal and rotational Killing fields. The effective potential is given by 



where k = —1,0,1 for timelike, null, and spacelike geodesies, respectively. We are most 
interested in the radial turning point of the geodesic, denoted by r^^^, which corresponds to 
the largest root of Vcs- In particular, we wish to ask, given generic /(r) and h{r) describing 
our asymptotically AdS spacetime, how small can we make rmin by adjusting k, E, and L? 
It is clear from the form of the potential that for fixed E and L, null geodesies have higher 
effective potential (outside the horizon) than spacelike ones, which in turn means that the 
corresponding Tmin will be larger for a null geodesic than for the corresponding spacelike one. 
This quick argument is further substantiated by actually determining where the turning 
point can occur. This is carried out in [209] ; here we will only summarize the main points: 
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• For causally trivial spacetimes without an event horizon, one can probe the entire space- 
time, both with spacelike and with null geodesies. In particular, for radial geodesies. 

However, here it is interesting to ask a more refined question: suppose we have access 
only to a certain region on the boundary, in particular we can only use geodesies with 
both endpoints in that region; out of this constrained set of geodesies, which one probes 
the deepest? While the answer depends on the details of the geometry, in AdS we can 
explicitly confirm that for a fixed angular span of the endpoints, rmin is minimized for 
E = (which necessarily implies spacelike geodesic) and grows monotonically with E"^. 
In fact, null geodesies probe the bulk least, as in the limit of E ^ oo the spacelike 
geodesic approaches a null one. 

• For spacetimes with an event horizon, null geodesies can probe only down to a finite 
distance outside the horizon, whereas specific spacelike geodesies (namely with E = 
and L approaching the horizon size from above) can probe arbitrarily close to the 
horizon. However, no geodesic with both endpoints anchored to the boundary can 
probe past the horizon. 

It is important to keep in mind that the above statements pertain to static spacetimes of the 
form (7.1), and do not necessarily hold for rapidly evolving spacetimes. In particular, if one 
collapses a black hole sufficiently quickly, spacelike geodesies can penetrate inside the event 
horizon. This can be easily seen by the simple gedanken experiment of [190], whose main 
point was to argue that precisely due to its teleological nature, the event horizon cannot 
pose a limitation to how far such bulk probes can reach. The idea, sketched in Fig. 9, is 
simply the following: Suppose we collapse null shell with sufficiently high energy to create 
a very large black hole. Inside the shell, the spacetime is that of pure AdS, so the spacelike 
geodesies just to the past of the shell creation probe the full spatial slice of the bulk; yet the 
event horizon from the black hole formed by the shell extends though this pure AdS region, 
and therefore gets penetrated by these unsuspecting geodesies (such as g in Fig. 9). 

Having classified the behaviour of geodesic probes, let us now turn to higher-dimensional 
surfaces. The motivation for considering such surfaces is that one has access to observables 
in the field theory whose holographic avatars involve such surfaces. The fully story for these 
higher dimensional surfaces is rather more involved, so we will mostly draw lessons from cases 
with sufficiently high degree of symmetry. In particular, we will keep the bulk spacetime 
static and we will consider a bulk minimal surface at constant time anchored on two parallel 
{n — l)-dimensional planes or an [n — l)-sphere on the boundary, with the boundary metric 
that of (i-dimensional Minkowski spacetime. To warm up we'll start with the bulk spacetime 

More specifically, in all cases the geodesies can only probe to depth given by circular orbit; for null 
geodesies this location is fully specified by the geometry, whereas for spacelike geodesies, by adjusting L we 
can bring this orbit arbitrarily close to the horizon. 
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Fig. 9: Gedanken experiment of [190] demonstrating that spacelike geodesies g (or minimal surfaces) 
which are anchored on the AdS boundary can pass through an event horizon H of a dynamically 
evolving spacetime, here exemplified by a null shell s collapsing to form the black hole (shaded 
blue). The corresponding Penrose diagram is sketched in Fig. 8. 

being tliat of pure AdS in Poincare coordinates, and then move on to a more general class of 
asymptotically AdS bulk spacetimes. While some of the results mentioned below are well- 
known and have been derived previously (see e.g. [208, 171]), the following presentation is 
based on [209] which contains a more systematic analysis. 

Let us start with pure AdS^+i in Poincare coordinates, 

ds'^ = \ i-dt^ + dx^ + dy'^ + dz"^) , (7.3) 

where z is the bulk radial coordinate with z = corresponding to the AdS boundary. One of 
the simplest cases to consider is an n-dimensional minimal surface anchored on the boundary 
of a "strip" of width Ax, with infinite extent in the other n—1 directions. The cross section 
z{x) of the corresponding minimal surface is obtained by minimizing its area. Since the 
area depends on n, the cross-section z{x) likewise depends on n, despite the translational 
invariance in n — 1 directions. Although one can easily obtain the cross section analytically, 
here we content ourselves with indicating its behaviour graphically. The left panel of Fig. 10 
shows the cross section of an n-dimensional minimal surface for n = 1,2, ... ,8, with the 
outermost curve representing n = 1, i.e. a spacelike geodesic, which has semi-circular shape. 
We see that as we increase the dimensionality, for a fixed depth z^, to which such a surface 
reaches, the strip width Ax decreases. Conversely, keeping Ax fixed, higher dimensional 
surfaces reach deeper into the bulk. This is summarized in the right panel of Fig. 10 where 
the ratio of depth probed to boundary extent Ax is plotted as a function of n; we see that 
it increases approximately linearly. 

The actual relation is somewhat more complicated, given by = ^ j^^tIttt- 



75 



zJAx 




Fig. 10: Left: cross-section of n-dimensional minimal surfaces in Poincare AdS^+i, with varying 
dimensionality n = 1,2, ...,8: the outermost (red) curve corresponds to n = 1 while the 
innermost (purple) curve to n = 8. Note that d (as long as it's high enough to accommodate 
the surface) does not enter. Right: corresponding ratio of maximal bulk radial extent to 
its boundary size Ax. 



One might conclude from tliis observation that higher-dimensional surfaces are therefore 
better probes of the bulk geometry. This will indeed partially motivate §7.3; but the indi- 
cation given above is a bit too glib. In particular, we have only considered the extent in the 
X direction and ignored the fact that the higher dimensional surfaces utilize extra directions 
of infinite extent. A slightly better comparison would therefore involve a boundary region 
of finite extent in all directions. The most natural such region is a (rf — l)-ball (with some 
specified radius Rq) in R*^"^'^. Let us therefore compare minimal surfaces anchored on S"^^ 
of radius Rq, for different values of n. Although the equation of motion for such a surface 
depends on n, this equation is satisfied by a very simple solution^^ z{f) = ^/ R^ — f'^ for all 
n - that is, the minimal surface is simply an n-hemisphere, whose bulk extent z^: = Rq is 
independent of n. However, it is clear that this can happen only due to a special cancellation 
in the bulk geometry; for general geometries the minimal surface profile z{r) would depend 
on n. The interesting question, then, becomes, for a fixed boundary region specified by Rq, 
how does the depth probed z=k vary with n? 

This question can be studied numerically for various spacetimes. Perhaps not surprisingly, 
it turns out that for "sensible" spacetimes, such as the planar Schwarzschild-AdS black 
hole, similar result as in Fig. 10 holds: the higher- dimensional surfaces again probe deeper. 
However, one can easily write a bulk spacetime metric violating energy conditions for which 
the opposite effect takes place: lower-dimensional surfaces probe deeper. Nevertheless, since 
such cases are unphysical, we will take the black hole the prototypical one. The 

basic rule of thumb is then very simple and consistent with naive expectations based on the 

Here f denotes the boundary radial variable, such that the corresponding boundary metric takes the 
form ds^ = [—dt^ + df^ + d^^_i + dy^ + dz^] , and by (n — l)-spherical symmetry, the minimal surface 
will be described by a single function z(f). 



76 



UV/IR duality: If one has access to only a certain region A on the boundary, and therefore 
is constrained to probes contained entirely within that region, then one can probe maximal 
amount of the bulk geometry by utilizing the full region A. In particular, if ^ is d — 1 
dimensional, and thus covers some open set in the boundary space, then one should use 
probes dual to the full d — 1-dimensional surfaces in the bulk for probing the bulk geometry. 

Finally, ere ending this section, let us revisit the issue of penetrating the horizon of a 
static bulk spacetime in the context of these higher-dimensional minimal surfaces. We have 
seen thtat geodesies with both endpoints anchored on the boundary cannot reach past a 
horizon, based on quite general arguments independent of the details of the bulk geometry. 
But we have also seen that higher-dimensional surfaces can probe deeper then geodesies; so 
one might naturally wonder whether such surfaces could in fact be used to probe past the 
horizon. The answer turns out to be negative, and similarly robust as for the geodesic case: 
no minimal surface anchored on the boundary can reach past the horizon. Therefore, in this 
sense, if we have access to the full boundary, for static spacetimes, spacelike geodesies are as 
good probes as higher-dimensional extremal surfaces, since in each case, the horizon is the 
limit to how deep these types of probes reach. Although we have so far focused on static 
geometry, shortly in §7.3 we will see an example of the behaviour of extremal surfaces for 
dynamically evolving geometry describing a black hole collapse. 

7.3 Entanglement entropy as probe of bulk dynamics 

Motivated by the discussion in §7.2, we are naturally led to consider spacelike extremal 
hypersurfaces anchored on the boundary with highest dimensionality. Ryu and Takayanagi 
[207, 208] proposed that in static geometries, the area of a bulk minimal co-dimension two 
spacelike surface (at a constant time) anchored on the boundary dA or certain region A 
on the boundary of AdS gives the entanglement entropy of the region A in the dual CFT 
state.^^ This was generalized in [171] to general time-dependent configurations, where the 
entanglement entropy is given by the area of an extremal surface, which is in fact related 
to light-sheet constructions of the covariant entropy bound [212] in the bulk spacetime. It 
is gratifying to note that such class of special highest-dimesionality surfaces are related to 
an appropriately special notion on the dual side: entanglement entropy is an important 
concept in field theory systems, with apphcations to condensed matter systems, quantum 
information, etc.. Further, given its non-extensive nature i.e., area scaling [213, 214], one is 
tempted, in the context of holography, to think of it as providing a measure for the effective 
degrees of freedom associated with a given region. This is indeed supported by the proposal 

Justification of this remarkable statement was subsequently provided in [210] using replica trick to 
evaluate entanglement Renyi entropies; although an error in this derivation was recently pointed out by 
[211], the correction nevertheless gives further evidence for the correspondence. 
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of [207, 208] and its covariant extension [171].^'^ 

To be slightly more specific, let us recall the definition of entanglement entropy and the 
holographic prescription for computing it. Consider a quantum field theory on dAi = A/'x 
and focus on a region A C Af. Given a density matrix p (or a pure state) on Af we define the 
reduced density matrix p_4 = Tr_^\_4(p). The entanglement entropy is then simply the von 
Neumann entropy associated with px- 5*^ = — Tr (p_4 log pX). The covariant holographic 
prescription [171] for computing this is to consider the problem of finding a co-dimension 
two extremal surface in the bulk geometry M. (dual to the state in question), which 
ends on the boundary of the chosen region A at the spacetime boundary dA4. S_a is then 
given by the area of Wa in Planck units i.e., S_a = In the situation where the 

bulk geometry is static (and hence the dual state is time-invariant) the problem reduces to 
the simpler one of finding minimal area surfaces [207, 208]. An excellent summary of the 
developments in understanding this holographic entanglement entropy can be found in [217]. 

Let us now apply this prescription to the interesting case of black hole collapse. In [171] 
a setting rather similar to the ones discussed earlier in §6.1.2 was considered to ascertain the 
time evolution of entanglement entropy. In particular, these authors considered the Vaidya 
solution in global AdS, 



with m{v) = tanhf smoothly interpolating between (corresponding to pure AdS) and 1 
(corresponding to Schwarzschild-AdS), in 3 bulk dimensions. The corresponding extremal 
surfaces (which in 2-1-1 dimensions are simply spacelike geodesies) are plotted in Fig. 11. 
At early times, before a black hole forms, these minimal surfaces probe the entire bulk, 
r G (0, oo). However, as the black hole forms and grows, these surfaces are repelled by the 
horizon, so that the region of spacetime they are sensitive to is only r G (r+(f ), oo). 

Having obtained the extremal surfaces for the Vaidya-AdS geometry, we can compute 
the entanglement entropy of the region A using the area of these surfaces. In the case 
of Vaidya- AdSs spacetimes one can in fact carry out this computation explicitly, since the 
extremal surfaces are simply spacelike geodesies in the geometry. The entanglement entropy 
is then just the regularized length Lreg of these geodesies. In fact, the calculation can be 
done analytically if the variation of the mass is slow, analogously to the discussion in §5.3.2. 
The numerical evaluation of the time-dependent entanglement entropy in this case and its 
comparison to the analytic approximation is shown in Fig. 12. For matter fields in the bulk 
satisfying the null energy conditions, and for a fixed boundary region, one can can show 
that the entanglement entropy increases (at least in the adiabatic approximation). Hence 

^^Further evidence of the entanglement entropy providing a measure of the active degrees of freedom is 
provided in [215] in the context of induced gravity (brane-world) models. See also [216] for recent attempts 
to reverse-engineer spacetime geometry from knowledge of entanglement entropy. 




(7.4) 



78 



vO = 0.1 vO = 0.5 vO = 1 




Fig. 11: Spatial projection of extremal surfaces in Vaidya-AdS (collapsing black hole) geometry, 
calculated at various times in the black hole formation process. The outer (black) circles 
correspond to the AdS boundary at time given by the vO indicated; the inner (red) circles 
correspond to the event horizon cross-section at corresponding time. The remaining curves 
correspond to the extremal surfaces anchored on the boundary, colour-coded by the size of 
the boundary region they encompass. At early times vO before the black hole forms, these 
extremal surfaces coincide with minimal surfaces in pure AdS. We can see that as the black 
hole grows, these surfaces are 'repelled' to remain outside the horizon. 




Fig. 12: Regularised proper length Lreg, or equivalently, the entanglement entropy S_a, as a function of 
the boundary Voo, for the particular region (an arc of fixed length along the spatial direction) 
in the Vaidya spacetime (red dots) and the corresponding analytic prediction in BTZ (black 
curve). 
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assuming that the matter undergoing collapse to form the black hole satisfies the null energy 
condition, it is evident from the adiabatic approximation that the entanglement entropy 
Sj[{v) oc Lregiy) iucreases in the process of a gravitational collapse. 

As suggested by our discussion above, it is clear geometrically that the entanglement 
entropy being described by an extremal surface should be sensitive to the causal structure of 
the spacetime. This is borne out explicitly in the computations of [171] which provide a strong 
coupling counterpart to the field theoretic results on temporal evolution of entanglement 
entropy [218]. 

In summary, we have seen that there exists a rich variety of bulk probes, with different 
merits. The bulk-cone singularities discussed in §7.1 probe smallest region of the bulk, but are 
quite robust and presumably easiest to extract as CFT 'data'. Spacelike geodesies or higher- 
dimensional minimal surfaces probe deeper, and as we've exphcitly seen in §7.2, can be used 
even when the boundary lives on a non-compact flat spacetime. Moreover, unlike bulk-cone 
singularities, these probes are all sensitive to the conformal factor of the bulk spacetime, 
but on the other hand, for static geometries, any surface at constant time is completely 
insensitive to the time component of the metric. Nevertheless, by considering finite-energy 
spacelike geodesies, or extremal surfaces which are not localized to constant time, we can 
probe the temporal aspect of the geometry as well. Although part of our discussion focused 
on static geometries because these allowed us to make more general statements and were 
easier to analyze, the most important point of this section in the present context is that 
these probes can tell us much about physics far out of equilibrium. 

8 Conclusion 

In this Review we have focussed on various aspects of non-equilibrium physics that can be 
accessed via holographic methods. The intriguing aspect of studying dynamical phenomena 
in AdS/CFT, apart from the intrinsic interest in developing computational techniques to 
facilitate investigation of strong coupling physics, is the natural mapping of these dynamics 
into the time dependent phenomena in a classical bulk spacetime. 

We have in the course of this Review touched upon various developments in the field, 
ranging from linear response regime, to the hydrodynamic regime, and finally to the fully non- 
linear regime. While aspects of the physics in various situations are by now well understood, 
it is clear that there are still open issues. In the following we will discuss some of the 
interesting open questions which could provide a fruitful avenue for further progress. 

In §4.3 we have mentioned some of the seminal attempts to model probes moving through 
a thermal medium, which is of relevance to the QGP phenomenology. The thermal medium 
is provided by a black hole geometry in the gravitational dual, while a probe in the form 
of an external quark is modelled by a fundamental open string (ending on the quark) in 
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this background. The motion of the quark determines the trajectory of the string's world- 
sheet, which in turn backreacts on the bulk metric to induce a non-trivial stress tensor on 
the boundary. The latter describes not only the bare quark, but also the cloud of partons 
around it, as well as the radiation produced by the quark's motion. It is interesting that these 
apparently distinct manifestations of the quark's motion are unified by a single fundamental 
quantity, namely the string, in the bulk. Conversely, given the bulk string, one might want 
to ask, which part of this string describes the quark and its surrounding halo, and which 
part the radiation?^*^ In this context, it has been argued in [102, 105, 106] building on the 
formalism developed in [219] that it is indeed possible to separate the total four- momentum 
of the string into a part intrinsic to the boundary quark, and another that corresponds to 
the radiation. The results of these computations seem to differ from the commonly assumed 
picture quoted in footnote 58, which has been advocated earlier in [220, 103, 221].^^ 

The fluid/gravity correspondence was reviewed in §5. From the geometrical standpoint, 
several intriguing issues deserve further investigation. Note that the fluid/gravity correspon- 
dence establishes a one-to-one map between long-wavelength gravitational solutions and 
solutions of fluid dynamics. Naively, one might think of this as merely a low-energy effective 
description of the AdS/CFT correspondence. However is it not true that any gravitational 
solution in AdS admits a fluid description. In particular, it would be useful to understand 
the role of non-long-wavelength bulk classical or semiclassical solutions, such as small AdS 
black holes. Likewise, there is room for more detailed bulk analysis, such as probing the 
allowed horizon topology and dynamics in more general situations, the nature of curvature 
singularity, or Cosmic Censorship. 

More generally, the fluid/gravity correspondence offers new insight into the black hole 
membrane paradigm, [222]. As described earlier, conventional membrane paradigm (cf., 
[135]) provides a simple picture of black hole dynamics in terms of classical physics of fluid 
living on a "membrane" (or stretched horizon) just outside the event horizon. Taking a more 
general view of trying to encode the black hole dynamics by fluid dynamics localized on a 
membrane in the spacetime, the immediate natural question is: where should such a mem- 
brane live? Perhaps the most obvious candidate is the event horizon; but this is problematic 
due to its null nature, and more importantly, because it is globally defined so we can't fix 
its position without knowing the full future evolution of the spacetime. Alternately, several 
(quasi)local notions of a black hole have been proposed, such as the so-called dynamical 
horizon, which however are spacelike surfaces inside the event horizon, and therefore do not 
admit the standard notion of evolution. A more popular suggestion is the stretched horizon, 
which is the formulation given by the membrane paradigm [135]. However, there likewise 
remain ambiguities in localizing stretched horizon. Within the fluid/gravity correspondence, 

^® It is in fact commonly assumed that the part inside the world-sheet horizon corresponds to the radiation, 
since this portion can no longer causally influence the quark, whereas the part above the horizon contributes 
to the halo. 

^^We thank Alberto Guijosa for alerting us to these references. 
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the full spacetime dynamics is mapped to the dynamics of a conformal fluid, which albeit 
reminiscent of the membrane paradigm, has one important twist: the membrane lives on the 
boundary of the spacetime (which is unambiguously defined and admits a fluid description 
with well-defined dynamics), and gives a perfect mirror of the bulk physics. This "membrane 
at the end of the universe" picture is a natural consequence of the holographic nature of the 
fluid/gravity correspondence. Recently this view has been countered in [223], where the 
authors argue that the projection of Einstein's equations onto any null hypersurface leads in 
a long-wavelength approximation to hydrodynamical behaviour. While this is interesting, it 
does not address the objections given above and more importantly it seems to us somewhat 
unnatural that a feature common to all null hypersurfaces has something fundamental to do 
specifically with black hole horizons.^" 

As mentioned in §1, one of the outstanding problems in fluid dynamics is understanding 
turbulence. We might therefore ask how relevant is the fluid/gravity correspondence (ap- 
propriately generalized) to tackling such a problem. We argue that the correspondence is 
indeed relevant and worth vigorously pursuing further (see [148, 158, 223, 225] for thoughts 
on these issues). The important point is that the regime where our derivative expansion is 
valid corresponds to large Reynolds number (where the viscous terms are small relative to 
the leading order terms), so the phenomenology of turbulence should be directly relevant 
to the study of near-equilibrium AdS black hole dynamics. At least for the widely stud- 
ied case of non-relativistic fluids described by the Navier-Stokes equations, turbulence has 
many striking qualitative features, including the tendency to dynamically break symmetries 
as Reynolds number is increased, the sharp onset of turbulence at critical Reynolds num- 
bers, and an "energy cascade" in fully developed turbulence in which energy is transferred 
in a predictable way between eddies at different scales. Less is known about turbulence for 
the microscopically relativistic fluids relevant in our context, but it would be fascinating 
to understand the gravitational interpretation for those features that do generalize to the 
relativistic case. 

Finally, one intriguing observation, made in §6.1.1, is that in certain cases of large dispar- 
ity between apparent and event horizon, such as in the conformal soliton setting analyzed in 
[170], the CFT dual entropy seems to be given by area of the apparent horizon rather than 
the event horizon as conventionally expected. Since this striking observation is rather sur- 
prising and, if true in general, would have far-reaching implications, let us end our discussion 
by probing this issue a bit deeper. In particular, can it be true that the entropy is associated 
with the area of the apparent horizon? The fact that CFT entropy is not determined by 
the area of the event horizon in general is perhaps not so surprising, given that the event 

^°More recently [224] argue that a Wilson renormalization group perspective allows one to relate the 
'membrane at the end of the universe' picture of the fluid-gravity correspondence with the conventional 
membrane paradigm story. In particular, the authors argue that as the location of the stretched horizon is 
brought in from the boundary to the actual horizon, one in fact sees an RG flow from a relativistic conformal 
fluid to a non-relativistic incompressible fluid. 
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horizon is a global construct - its precise location behaves teleologically: it depends on the 
spacetime arbitrarily far in the future - whereas CFT entropy at a given time is measured 
independently of what happens at a later time.^^ On the other hand, there are nevertheless 
serious problems with associating the area of an apparent horizon, or any other quasilocal^^ 
horizon, to the dual CFT state entropy. One difficulty is that the entropy is expected to 
be smoothly varying in time, whereas the area of a given apparent horizon can jump dis- 
continuously. A more problematic issue, however, is that the location of apparent horizons 
depends on a choice of spacetime foliation, whereas the CFT state at a given time carries 
no such information. This seems to imply that our bulk prescription has more freedom or 
ambiguity in defining the entropy than that afforded by the boundary theory. 

One possible resolution is that there is a preferred foliation of the bulk, such as a zero- 
mean-curvature slicing, on which one is supposed to evaluate the area. However, we have no 
compelling justification for this option. A simpler resolution to this puzzle would be that in 
the regime where the concept of entropy is meaningful, the horizon has to be evolving slowly 
enough that there is negligible difference between the areas of all slices of horizon which end 
on the same boundary time-slice. This is essentially the same picture as that advocated in 
[22, 149], except that here we use it for apparent horizon rather than the event horizon. In 
effect the field theory on the boundary should achieve local equilibrium in order for entropy 
to constitute a meaningful observable. 

However, intriguingly even this proposition cannot be the full story. If there is no event 
horizon at early times, such as in the case of a black hole formed from a collapsed star, then 
we can find foliations without apparent horizon for arbitrarily late times. A simple example 
of such a foliation is sketched in Fig. 13. Suppose we collapse a black hole, and wait for very 
long time after the black hole has settled down, say time tiate indicated in Fig. 13. A 'natural' 
foliation (sketched in the left panel) would give the expected results: by the time tiate, the 
apparent horizon area has grown and settled down to a constant which coincides with the 
event horizon area. However, we can devise a foliation (sketched in the right panel) which, 
even at tiate, admits no trapped surfaces and therefore has zero apparent horizon area. The 
corresponding spacelike surface 'dips just below' the event horizon formation; and since the 

Although some might argue that there is httle distinction between the two cases, since with certain 
analyticity assumptions on the spacetime the evolution of the entire spacetime is determined by the initial 
state (and the boundary conditions), we do not take this viewpoint: there is a difference between simply 
measuring certain properties of a state at a given time, and evolving to that state from much earlier initial 
time - the latter requires far more detailed information. 

For a good review of various notions of quasilocal horizons, see [226, 227]. Technically speaking, an 
apparent horizon is defined as a co-dimension 2 surface, on a given leaf of foliation, corresponding to the 
outermost marginally trapped surface or the boundary of trapped points. In fact, the set of apparent horizons 
on all leaves of a given foliation need not even form a smooth co-dimension 1 surface in the full spacetime. 
In the present review, we will however stick to the more commonly used terminology in non-GR literature 
of apparent horizon signifying the full timclike co-dimension 1 surface. 



83 




Fig. 13: Possible spacelike foliations of a geometry a black hole collapse. The outer cylinder represents 
the AdS boundary; the bottom (orange) surface represents a collapsing star, the (red) vertical 
line the curvature singularity, the inner (purple) cylinder the event horizon. Two distinct sets 
of foliations are sketched. On left, the corresponding apparent horizon has settled down at 
late time tiatc (indicated in blue). On right, the spacelike surface (green) at tiatc admits no 
apparent horizon at all. 

event horizon itself is null, the requisite foliation surface can always be constructed so as to 
be spacelike. 

It would be interesting to understand the full implications of this observation, since it 
bears on the discussion in §1 regarding the relation between boundary time and bulk time: 
which bulk region is most naturally associated with a given boundary time, such as tiate? 
Since the notion of entropy depends on amount of coarse-graining, perhaps this is the natural 
feature which the choice of bulk foliation reflects. 
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